UNIQUENESS RESULTS FOR MATRIX- VALUED 
SCHRODINGER, JACOBI, AND DIRAC-TYPE OPERATORS 

FRITZ GESZTESY, ALEXANDER KISELEV, AND KONSTANTIN A. MAKAROV 

Abstract. Let g(z,x) denote the diagonal Green's matrix of a self-adjoint 
m X m matrix- valued Schrodinger operator H = —-^jlm + Q(x) in L 2 (R) m , 
m £ N. One of the principal results proven in this paper states that for a fixed 
xo £ M and all z £ C+, g(z,xo) and g'(z,xo) uniquely determine the matrix- 
valued m X m potential Q(x) for a.e. We also prove the following local 
version of this result. Let gj(z,x), j = 1,2 be the diagonal Green's matrices 
of the self-adjoint Schrodinger operators Hj = — -j^Im + Qj{&) m L 2 (R) m . 
Suppose that for fixed a > and xq £ M, g^z.xo) — g2(z, xo)\\ cmX m + 
||gi(z, xo) — g' 2 (z, xo)\\ C mxm = 0(e _2Im ( z ^ ) a ) for z inside a cone along 

| 2 | — >0O 

the imaginary axis with vertex zero and opening angle less than it/ 2, excluding 
the real axis. Then Qi(x) = Q 2 (x) for a.e. x £ [xq — a, xq + a]. 

Analogous results are proved for matrix-valued Jacobi and Dirac-type op- 
erators. 

This is a revised and updated version of a previously archived file. 



1. Introduction 

While various aspects of inverse spectral theory for scalar Schrodinger, Jacobi, 
and Dirac-type operators, and more generally, for 2x2 Hamiltonian systems, are 
well-understood by now, the corresponding theory for such operators and Hamil- 
tonian systems with m x m, m 6 N, matrix-valued coefficients is still in its infancy. 
A particular inverse spectral theory aspect we have in mind is that of determining 
isospectral sets (manifolds) of such systems. It may, perhaps, come as a surprise 
that determining the isospectral set of Hamiltonian systems with matrix- valued pe- 
riodic coefficients is a completely open problem. It appears to be no exaggeration to 
claim that (unless one considers trivial cases such as diagonal coefficient matrices, 
etc.) absolutely nothing seems to be known about the corresponding isospectral 
sets of periodic matrix-valued Schrodinger operators with the sole exception of its 
compactness. The same ignorance applies to Jacobi, Dirac, and more generally, to 
periodic 2m x 2m Hamiltonian systems with m > 2. While the present paper is 
not sufficiently ambitious to change this sorry state of affairs, we will take a mod- 
est step toward a closer investigation of inverse spectral problems and prove a few 
uniqueness theorems for such systems, that is, determine spectral data (interpreted 
in a very broad sense) that uniquely determine the matrix-valued coefficients in 
Schrodinger, Jacobi, and Dirac-type systems. It should be mentioned that these 
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types of problems are not just of interest in a spectral theoretic context, but due 
to their implications for other areas such as completely integrable systems (e.g., 
the nonabelian Korteweg-deVries, Toda lattice, and Ablowitz-Kaup-Newell-Segur 
(nonlinear Schrodinger) hierarchies), are also of considerable interest to a much 
larger audience. 

Before we continue along this line of thought, it seems appropriate to briefly 
mention some of the present day knowledge of 2m x 2m matrix- valued Hamiltonian 
systems and the Weyl-Titchmarsh and spectral theory associated with them. To 
save space we will often simultaneously discuss references on Schrodinger, Dirac, 
and general Hamiltonian systems together without differentiating between them, 
and occasionally single out the finite-difference (Jacobi) systems. Moreover, there 
exists a considerable amount of literature on the foundations of Hamiltonian sys- 
tems (and their special cases, such as Schrodinger, Jacobi, and Dirac-type systems) 
which necessarily forces us to be rather selective here. In particular, we focus 
primarily on the case m > 2 for the remainder of this introduction. 

The basic Weyl-Titchmarsh theory of regular Hamiltonian systems can be found 
in Atkinson's monograph |T^| ; Weyl-Titchmarsh theory of singular Hamiltonian 
systems and their basic spectral theory was developed by Hinton and Shaw and 
many others (see, e.g., 0Jf§-f7|, @, @, @, ||, jT^, fl8|, 

122 , |12S|1 , 132 , Ch. 9], |l40| and the references therein); the corresponding the- 
ory for Jacobi systems can be found in |Q, [13C] and the literature therein. 
Various aspects of direct spectral theory, including investigations of the nature of 
the spectrum involved, (regularized) trace formulas, etc., appeared in |lT 

a is 



§, pg, |§, BJ, M, Kg, [109|, [114|, hlM. General asymptotic 



expansions of Weyl-Titchmarsh matrices as the (complex) spectral parameter tends 
to infinity under optimal regularity assumptions on the coefficients are of relatively 

p37|). The inverse 



(see also [126 



recent origin and can be found in 
scattering formalism for continuous Hamiltonian systems has been studied by a va- 
riety of authors and we refer, for instance, to§,J§,§,@,0>IIl' ©> 



111 , |138[], [ 145 1 . General inverse spectral theory, the existence of transformation 



operators, etc., is discussed in (62), @, @, gOj, p21j , pjfl , p5| , |2§, p9fl , 
131 , [ 132 1 , and the references therein. Inverse monodromy problems for canoni- 
cal systems received a lot of attention recently. The interested reader is referred 



to 0, @, Hi, |^-[jK)2|], p7| , p3^ | and the extensive literature cited therein. 
The corresponding inverse spectral and scattering theory for matrix-valued finite 
difference systems and its intimate connection to matrix-valued ortho gon al poly- 
nomials is treated in (|, §|, @, @, 



l-m 



97 



|, |130(| . More 

specific inverse spectral problems, such as compactness of the isospectral set of 
periodic Schrodinger operators [Q, special isospectral matrix- valued Schrodinger 
operators, and Borg-type uniqueness theorems (for periodic coefficients as well as 
eigenvalue problems on compact in t erva ls) wer e recently studied in pH , J32|, |37j, 

3j§, JH, |7|, j7|, ||, fiol , fiol , H|, |||. Moreover, direct spectral theory in 
the particular case of periodic Hamiltonian systems (i.e., Floquet theory and alike) 
has been studied in j§, |tJ, f§, @, g|, g, g, ||, g, gT^, (u|-@, 
with many more pertinent references to be found therein. 

Apart from Floquet theoretic applications in connection with Hamiltonian sys- 
tems already briefly touched upon, we also need to mention applications to random 
Schrodinger operators associated with strips as discussed, for instance, in [ p9| , p2| , 

87 1, Q, and especially to nonabelian completely integrable systems. Since the 
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literature associated with the latter topic is of enormous proportions, we can only 
refer to a few pertinent publications, such as, O, §(!], {HI, j3pj, @-@], |1, 



103 , 107 , p8[, 112 , p3|, [126 -|128|, [13C]. The interested reader will find a 



wealth of additional material in these references. 

Finally, we turn to the principal subject of this paper, that is, uniqueness-type 
theorems for Schrodinger, Jacobi, and Dirac-type operators. We have already men- 
tioned a few uniqueness results in connection with Borg-type theorems for such sys- 
tems. Additional uniqueness-type results in terms of matrix-valued Wronskians and 
transformation operators can be found in some papers by Leibenzon [po| , plfl and 



Malamud ]100| |, 1 101 ] . The uniqueness theorems proven in this paper are directly 
formulated in terms of diagonal Green's matrices g(z,xo) and their x-derivatives 
g'(z,xo) at some fixed xq G R, for matrix- valued Schrodinger and Dirac-type op- 
erators on R and similarly in terms of diagonal Green's matrices g(z,ko,ko) and 
off-diagonal Green's matrices G(z,fc ,fco + 1), etc., for some fixed k G Z, in the 
context of matrix-valued Jacobi operators on Z. Moreover, we prove certain local- 
ized versions of these uniqueness theorems for exponentially close diagonal Green's 
matrices with respect to z as \z\ — * oo in the Schrodinger and Dirac-type context 
and analogous theorems in the case of Jacobi operators whose (diagonal and certain 
off-diagonal) Green's matrices differ by an inverse power of z as \z\ — - > oo. To be 
a bit more specific, we briefly describe some of our principal results in the case 
of self-adjoint Schrodinger operators H in L 2 (R) m , m G N associated with m x m 
matrix-valued differential expressions — -£pI m + Q(x), with Q = Q* G Lj L oc (R) mxm . 

Let C + be the open complex upper half-plane. Denoting by g(z,x), z G C\R, 
x G R the diagonal Green's matrix associated with H (i.e., the integral kernel of 
the resolvent of H on the diagonal, g(z,x) = (H — z)^ 1 (x,x)) and by g'(z,x) its 



x-derivative, we will prove the following result (cf. Theorem 3.7). 

Theorem 1.1. Fix xq G R. Then g(z,xo) and g'(z,xo) for all z G C+ uniquely 
determine the matrix-valued m x m potential Q(x) for a.e. x G R. 

Moreover, let Hj, j = 1,2 be two Schrodinger operators associated with the 
m x m matrix- valued potentials Qj{x) and gj(z,x) the corresponding diagonal 
Green's matrices of Hj. Then we prove the following local version of Theorem [0] 



(cf. Theorem 3. 



Theorem 1.2. Let a > and x G R. // 

\\gi(z,x ) - g 2 {z,xo)\\c™xm + ||^(z,a;o) - g' 2 (z, x )\\ C m.xm 

= 0(e- 2Im (* 1/2 > a ) (1.1) 

\z\ — >oc 

along a ray with < arg(z) < tt, then 

Qi(x) — Qzix) for a.e. x G [xo — a, xq + a]. (1.2) 

Our results were of course inspired by analogous ones in the special scalar context 
m = 1, but also by a very interesting uniqueness theorem proven by Berezanskii 
in 1953 ]Ti| , [ fl8| in the context of multi-dimensional Schrodinger operators H. 
Berezanskii's point of departure is a bit different from ours. He considered the 
spectral kernel #(A, x, y) of H (i.e., the integral kernel of the spectral projection 
Eh(X) of H , $(A, x, y) — Ejj(X, x, y)) and appropriate normal derivatives of it with 
respect to x and y across an arbitrarily small piece of a smooth surface in R d , 
d = 2, 3. There are additional hypotheses in Berezanskii's work which need not 
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be discussed here. In a certain sense we replaced the spectral projection of H by 
its resolvent and (in our one-dimensional context) the arbitrarily small piece of a 
surface by the point xq and the corresponding normal derivative by d/dx. 

In Section [2] we review the basic Weyl-Titchmarsh theory for matrix-valued 
Schrodinger and Dirac-type operators as needed in our principal Section [| The 
latter contains our uniqueness results for Schrodinger and Dirac-type operators. 
The final Section ^ then treats Weyl-Titchmarsh theory and uniqueness theorems 
for matrix-valued Jacobi operators. 

2. Matrix- Valued Schrodinger and Dirac-Type Operators 

In this section we briefly recall the Weyl-Titchmarsh theory for matrix-valued 
Schrodinger and Dirac-type operators. In order to treat both cases in parallel, we 
use the fact that both are special cases of Hamiltonian systems and hence develop 
the theory from that point of view. Throughout this paper all matrices will be 
considered over the field of complex numbers C and the corresponding linear space 
of k x I matrices will be denoted by C kxe . 

The basic assumption for Sections || and ||| of this paper will be the following. 

Hypothesis 2.1. Fix to G N and define the 2m x 2m matrix 
We consider two principal cases: Either 

(i) suppose Q = Q* E Ll oc (M.) mxm and introduce the 2m x 2to matrices 

0). m -{-«*> I). <-) 

or 

(ii) suppose 

A = I 2m , B = B* e Ll c (R) 2mx2m . (2.3) 

Given Hypothesis |2.1| we consider the Hamiltonian system 

M'(z, x) = (zA + B(x))&(z, x) for a.e. ieR, (2.4) 

where z G C plays the role of a spectral parameter and $f(z, x) is assumed to satisfy 

9{z, ■) G ACi oc (R) 2mx2m . (2.5) 

Here and later on, I p denotes the identity matrix in C pxp for p G N, M* the 
adjoint (i.e., complex conjugate transpose), M f the transpose of the matrix M, 
and ACioc(R) denotes the set of locally absolutely continuous functions on R. At 



times it will be convenient to consider a 2to x r solution matrix of (2.4), with 
r = 1, . . . , 2to, whose entries will then be assumed to lie in ACi oc (R)- 



One verifies that Hypothesis 2.1 (i) governs the case of matrix- valued Schrodinger 
operators. In fact, inserting 

*M=fe!'!i (2-6) 



into equation (2.4) then yields 

-ijj"(z, x) + Q(x)ipi(z, x) = zipi(z, x), (2.7) 
ip 2 (z,x) = ip[(z,x). (2.8) 
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Here it is assumed that 



^{z, •) e AC loc (Kr 



J = 1,2. 



(2.9) 



In the case of Hypothesis 2.1 (ii), (2.4) represents a Dirac-type system. 

In connection with our uniqueness result for Dirac-type operators in Section |^ 
we will make the following additional assumption due to the freedom of certain 
gauge transformations in the Dirac context, which leaves the corresponding spectral 
matrix, but not the potential B(x) invariant, as discussed by Gasymov ]49| ] and 
Gasymov and Levitan pi . 



Hypothesis 2.2. Suppose Hypothesis 2.1 and in the case of Dirac-type opera- 
tors {i.e., assuming (2.3)) suppose that B is locally essentially bounded, B G 
L^ c (M.) 2mx2m , and that B(x) is of the special form 

B( X) = ( B ^f\ v?V). ( 2j °) 

\B 1 ^{x) ~Bi A (x)J 
with B h i(x) = B h i(x)*, B h2 (x) = B ia {x)\ j = 1,2. 

Equation ( [2.1C ) represents a typical normal form inspired by the requirement 

JB(x) + B(x)J = for a.e. (2.11) 

Various possible normal forms for Dirac-type operators are discussed in p2[ , 
M, pi), pi, and in the monographs @, Ch. 9], M Ch. 7], [EM p. 193- 



195]. The continuity assumption on B in Hypothesis 2.2 is further discussed in the 



paragraph following Theorem 2.12 



Next we briefly turn to Weyl-Titchmarsh theory associated with (2.4) and recall 
some of the results developed by Hinton and Shaw in a series of papers devoted 
to spectral theory of (singular) Hamiltonian systems (see also Q, Q). 

While they discuss ( |2.4| ) under much more general hypotheses on A(x) and B(x), we 
confine ourselves here to the special c ases of matrix- valued Schrodinger and Dirac- 
type systems governed by Hypothesis 2A . The following facts on Weyl-Titchmarsh 
theory are taken from |36| and |37f] and hence we omit the corresponding proofs. 

Let ^(z, x, Xo) be a normalized fundamental system of solutions of (2.4) at some 
.To G K, that is, ^(z,x,xq) satisfies 

JV'(z,x) = {zA + B(x))V(z,x), zeC 

for a.e. and 



(2.12) 
(2.13) 



y(z,x ,x ) 
Moreover, we partition ^{z, x, xo) as 
$>(z,x,xo) = (^j.kfo^aJo))^^ 

'di(z,x,x ) 4>i(z,x,x )^ 

9 2 (z,x,x ) (f> 2 {z,x,x ) / 

where 6j(z,x,xo) and <pj(z, x, xq) for j = 1,2 are m x m matrices, entire with 
respect to z G C, and normalized according to (2.13), that is, 

9i(z,x ,x ) = <p 2 (z,x ,x ) = I m , <px(z,x ,x ) = 9 2 {z,x ,x ) = 0. (2.15) 



= (Q(z,x,x ) 3>(z,x,x )) = \\ 



(2.14) 



Next, let & G C" 



J 



1,2, introduce the matrix (3 — {f3\ (3 2 ) G 



l , and 



assume that 

rank(/3)=m, [3/3* = I m , and either lm(/3 2 0l) < or Im03 2 /?D > 0. (2.16) 
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One can then prove the following result. 

Lemma 2.3 (see, e.g., p^ ). Assume Hypothesis \2.\, le t Q(z, x, Xq) and <&{z, x, Xq) 
be defined as in (|2.14 ) ; and suppose (3 satisfies ( 2.16 ). Then, for c G R\{xo}, 
f3<&(z,c,Xo) is singular if and only if z is an eigenvalue for the regular boundary 
value problem given by ( |2.4| ), ( |2.5[ ) together with the separated boundary conditions 

(I m 0)W(z,x o ) = 0, f3#(z,c)=0. (2.17) 



For the regular boundary value problem described in Lemma 2.3, the boundary 
conditions in ( p. 17 ) are self-adjoint whenever Im(/32/3i) = 0. 

Lemma |2.3| provides appropriate conditions for defining a certain meromorphic 
to x to matrix M(z, c, Xq, (3). 

Definition 2.4. Assume Hypothes is |2.1| and let Q(z, x, xo), and <fr(z,x,xo) be de- 
fined as in (2.14) with [3 satisfying (2.16). For c ^ xq, and (3<&(z, c, xo) nonsingular 
let 

M(z,c,x ,f3) = -[P®{z : c,x )]- x \pQ{z : c,x )\. (2.18) 
M(z, c, Xq, (3) is said to be the Weyl-Titchmarsh M -function for the regular bound- 



ary value problem described in Lemma 2.3 



The Weyl-Titchmarsh Af-function in (2.18) is an to x to matrix- valued function 
with meromorphic entries whose poles correspond to eigenvalues for the regular 
boundary value problem given by (2.4), (2.5), and (2.17). Moreover, if M G C" 1 *"^ 
and one defines 



U(z,x,xq) 



ux(z,x,x ) 

U 2 (z, X, X ) 



*(z, X, Xq) 



(2.19) 



J 



with Uj(z, x, Xo) G C m 
condition at x = c in (2.17) whenever M 



1,2, then U(z, x, Xq) will satisfy the boundary 
M (z, c, Xo, /3). Intimately connected 



with the matrices introduced in Definition 2.4 is the set of to x to complex matrices 
known as the Weyl disk. 

To describe this set, we first introduce the matrix- valued function E(M): With 
c 7^ xo, z G 
M e C mxm , let 



where 



a(s,t,z) = 



, and with U(z,c,Xo) defined by (2.12) in terms of a matrix 
E(M) = a(x Q ,c,z)U(z,c,xo)*(iJ)U(z,c,xo), (2.20) 
(s - i)Im(z) 



|( S -t)Im(z)|' 



a{s,t) = a(s,t,i), a(z) = a(l, 0, z), 



(2.21) 



with s, t G M, s ^ t. 



Definition 2.5. Assume Hypothesis 2.1 and fix xq G K, c G M\{xo}, and z G C\K. 
Then D(z, c, x ) will denote the collection of all M G C mxm for which E(M ) < 0, 
where E(M) is defined in ( 2.18 ). D(z,c,xq) is said to be a Wey/ dzsfc The set of 
M G C mxm for which E(M) — is said to be a Weyl circle (even when m > 1). 

This definition leads to a presentation that is a generalization of the description 
first given by Weyl [141]; a presentation which is geometric in nature, involves the 
contractive matrices V G C mxm , such that VV* < I m , and provides the justi- 
fication for the geometric terms of circle and disk and (see, e.g., |6S|] , [S^j , [0, 

H). 



UNIQUENESS FOR SCHRODINGER, JACOBI, AND DIRAC-TYPE OPERATORS 7 

The disk has also been characterized in terms of matrices which statisfy hypothe- 
sis (2.16) and which serve as boundary data for the regular boundary value problem 
described in Lemma |2~3| (see, e.g. Q). In particular, V(z,c,xq) denotes the 
collection of all M G C mxm obtained by the construction given in (2.18) where 
c ^ xo, z G C/R, where (3 is the m x m matrices defined in hypothesis (2.16) for 
which a(c, xo, z)Im(/32/3i ) > 0. 

We note that the Weyl circle corresponds to the regular boundary value problems 
in Lemma 2.3 with separated, self-adjoint boundary conditions. 

Lemma 2.6 (@, ||, @). Let M G C mxm , c ^ x , and z G C\M. TTien, 
= if and on/y ttere is a (3 G C mx2m satisfying (3J(3* = and = Z m 
smc/i i/iai 

= /3U(z,c,x ), (2.22) 

where U(z,c,xo) is defined in ( ^.12| ) m terms of M. With f3 so defined, 

M = -[/?$(*, c, xo)]- 1 ^ Q(z, c, x )], (2.23) 

that is, M = M(z, c, xo, 13). 

Next, we recall a fundamental property associated with matrices in T>(z, c, xo). 

Lemma 2.7 (see, e.g., |37j]). //Af G V(z,c,x ), then 

a(c,x ,z)Im(M) > 0. (2.24) 

Moreover, whenever (3 G C mx2m satisfies (3 J (3* = and = I m; 

M(z, c, x , 13) = M(z, c, x , (3)*. (2.25) 

For c > Xo, the function M(z, c, x , (3), defined by ( [2.18 ), and satisfying (2.24), 
is said to be a matrix-valued Herglotz function of rank m. Hence, for Im(/?2/3i) = 
0, poles of M(z, c, xo, 13), c > xq, are at most of first order, are real, and have 
nonpositive residues. Such functions admit representations of the form 



M(z,c,x ,P) =d+ zC 2 + / dn(X,c,x ,/3) 



1 



A 



A - z 1 + A 2 



C > Xo, 

(2.26) 



where (with (3 fixed) C*2 > and C\ are constant m x m self-adjoint matrices, and 
f2( A, c, Xq, f3) is a nondecreasing m x m matrix-valued function such that 



||dfi(A,c,x ,/3)|| Cm x m (l + A 2 )- 1 <oo, 



(2.27) 



S1((A, fi], c, xo, (3) = limlim — / dv a(c, xo, v + ie)lm (M(j/ + it, c, xo, (3)) . 

(2.28) 

For the self-adjoint boundary value problems in connection with Schrodinger, Dirac, 
and Jacobi operators discussed in this paper, C*2 = in ( 2.26| ) (and also later 
in (2.41) and ( 2.55| )) and f2(A, c, xo, f3) is then piecewise constant with jump dis- 
continuities at the eigenvalues of the problem. Analogous statements apply to 
—M(z, c, xo, (3) if c < xo. 
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We further note that the sets D(z, c, xq) are closed, convex, and are nested with 
respect to increasing or decreasing values of c (cf. Q, j70|, (jn|, EjJ, [113]), that 
is, 



D(z, C2,Xq) C D(z, c\,xq) for x < c% < c 2 or c 2 < c\ < x . 



(2.29) 



Hence, the intersection of this nested sequence, as c — ^ ±oo, is nonempty, closed 
and convex. We say that this intersection is a limiting set for the nested sequence. 



Definition 2.8. Assume Hypothesis |2.1| and let D±(z,xo) denote the closed, con- 
vex set in the space of to x m matrices which is the limit, as c — > ±oo, of the 
nested collection of sets D(z, c, xq) given in Definition 2.5. D±(z,xo) is said to be 
a limiting disk. Elements of D±(z,xq) are denoted by M±(z,xq) £ C mxm . 



In light of the containment described in (2.29), for c ^ xq and z £ 

D±(z,x Q ) C D(z,c,x ). (2.30) 



Given its characterization in the second paragraph following Definition 2.5, we see 
that D±(z, xq) consists of those matrices M{z, c, Xq,/3) where (3 satisfies ( 2.1 6j ) . In 
particular, 



M±(z,x ) = M(z,c,x ,(3) 
for an appropriate choice of /3. 



(2.31) 



When D±(z, xo) is a singleton matrix, the system (2.4) is said to be in t he limit 
point (l.p.) case at ±oo. When D±(z,xq) has nonempty interior, then (2.4) is said 
to be in the limit circle (I.e.) case at ±oo. Indeed, for the case to = 1, the limit 
circle case corresponds to D±(z,xq) being a disk in C. 

These apparent geometric properties for the disk correspond to analytic proper- 
ties for the solutions of the Hamiltonian system (2.4), (2.5). To recall this correspon- 
dence, we introduce the following spaces in which we assume that — oo < a < b < oo, 



L 2 A ((a,b)) = \cj ) :(a,b) 



dx(4>{x),A(t>{x))c^ < oo j, (2.32a) 

N(z, oo) = {<;6 £ L 2 A ((c, oo)) | Jfi = (zA + B)cf> a.e. on (c, oo)}, (2.32b) 
N(z, -oo) = {</>£ L^((-oo, c)) | Jef) = (zA + B)<f) a.e. on (-oo, c)}, (2.32c) 

for some c £ R and z £ C. (Here (cf), ip)c = Y^j=i fij^j denotes the standard scalar 
product in C", abbreviating x G C" by x — {xi, ■ ■ ■ > XnY ) Both dimensions of the 
spaces in ( |2.32b| ) and (2.32c), dimc(iV(z, oo)) and dimc(A r (z, — oo)), are constant 
for z £ C± = {C £ C ±Im(C) > 0} (see, e.g., @, Q). One then observes that 
the Hamiltonian system (2.4) is in the limit point case at ±oo whenever 

(2.33) 



(2.34) 



dim<c(-ZV(z, ±oo)) = m for all z £ 
and in the limit circle case at ±oo whenever 

dimc(iV(z, ±oo)) = 2to for all z £ C 



Next we recall the fact that the Dirac-type systems considered in this paper are 
always in the limit point case at ±oo. 



Lemma 2.9. Assume Hypothesis 2.1. Then the limit point case holds for Dirac- 



type systems [i.e., for A{x) = hm in (2.4)) at ±oo 
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Lemma 2.£ , under varying sets of assumptions on B(x), is well-known to experts 
in the field. For instance, in the case m = 1 and with B\p(x) — B2,i(x) this fact 
can be found in [13S]. For B £ C(R) 2mx2m and a more general constant matrix 



A, this result is proven in j92j (their proof, however, extends to the current B £ 
i 1 1 oc (R) case). More generally, multi-dimensional Dirac operators with L 1 2 oc (M™)- 



type coefficients (and additional conditions) can be found in |95 . A short proof 
of Lemma |2.9| can be found in p7| . After completion of this paper we became 
aware of a recent preprint by Lesch and Malamud fl93[ which provides a thorough 
study of self-adjointness questions for more general Hamiltonian systems than those 
considered in this paper. 

In either the limit point or limit circle cases, M±(z, Xq) £ 8D±(z,xq) is said 
to be a half-line Weyl-Titchmarsh matrix. (In the l.p. case one of course has 
D±(z,xq) = 8D±(z,xq).) Each such matrix is associated with the construc- 
tion of a self-adjoint operator acting on L\([xo, ±00)) (~l AC([xo, ±00)) for the 
Hamiltonian system (|2.4|), ( |2.5| ). However, for those intermediate cases where 
m < dimc(N(z, ±00)) < 2m, Hinton and Schneider have shown that not every 
element of the boundary dD±(z,xo) of D±(z,xq) is a half-line Weyl-Titchmarsh 
matrix, and have characterized those elements of the boundary that arc (cf . p6| , 

• 

For later reference we summarize the principal results on M±(z,xq) in the fol- 
lowing theorem. 



Theorem 2.10 (|61[, |68|, |6J], |]72j, |82|). Assume Hypothesis |U| and suppose 
that z £ C\R, xq £ R. Then 

(i) ±M±(z, Xq) is an m x m matrix-valued Herglotz function of maximal rank. In 
particular, 

lm(±M ± (z,x )) > 0, zeC+, (2.35) 

M±(z,x ) = M±{z,x )*, (2.36) 

rank(M±(z, xo)) = m, (2.37) 

limM±(A + i£,xo) exists for a.e. X £ M. (2.38) 

Local singularities of ±M±(z,xq) and ^M±(z,xq)~ 1 are necessarily real and at 
most of first order in the sense that 

=Flim(ieM±(A + ie,x )) > 0, A £ M, (2.39) 

ej.0 



eio \M±(\ + ie,x ) J - ' y J 

(ii) ±M±(z, xo) admit the representation 

±M±(z,x ) = F±(x )+ /dr» ± (A,xo)((A-z)- 1 -A(l + A 2 )^ 1 ), (2.41) 

where 

F ± (x )=F ± (x )*, I ||dr! ± (A, a;o )|| Cm x m (l + A 2 )- 1 <oo (2.42) 



n±((A, ii] , x )= limlim- / dvlnv{±M±(v + ie,x )). (2.43) 

S10 elQ 7T J x+S 
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(iii) Define the 2m x m matrices 



9i(z,x,x ) <t>i(z,x,x ) 



then 



K 6 2 (z,x,x a ) <j} 2 (z,x,x )J \M±(z,xq) 

/-±oc 

lm(M±(z, xq)) = lm(z) / dx *&±(z, x, xq)*A^/±(z, x, xq) 



(2.44) 



(2.45) 



where A has been introduced in (2.2), (2.3). 

In order to describe the Green's matrix associated with J2.12 ), we assume the 
hypotheses of Theorem [2.10| and introduce 

K(z, x, x') = v E'=p(z, x, Xo)[M_(z, xq) — M + (z, Xo)] "9±(z, x', xq)* , 

x § x', z e C\R, (2.46) 

and 

M{z,x ) = (M jd ,(z,x )).., =12 

= [K{z,x ,Xo + Q)+K(z,x ,x -0)}/2, z G C\R, (2.47) 
M hl (z,x ) = [M_(z,x ) - M + (z,x Q )}-\ 

Mi, 2 (z, x ) = 2- 1 {M_(z, x Q ) - M+{z, xo)]- 1 ^^, x ) + M+(z, x )}, 
M 2tl (z,x ) = 2- 1 {M-(z,x ) +M+{z,x )}[M-(z,xo) - M+foaso)] -1 , 
M 2 ,2(z,x a ) = M±(z,xo)[M-(z,xo) - M + (z, xq)]' 1 M T (z, x )- 



Next let <j> £ L\(M.) and consider 



Jip'(z, x) = (zA + B(x))ijj(z, x) + A<j)(x), z G 



(2.48) 



for a.e. x G M. Then ( p.48[ ) has a unique solution ip(z, ■ ) <E L 2 A (R) n ACi oc ( 
given by @, @ 



ip(z,x) 



dx' K(z,x,x')A(j)(x'). 



In the following we associate an operator T in L 2 (M.) m , respectively, L 2 ( 
the Hamiltonian system (2.4) denoted by 



T = 



H if A = 



Im 



/ 

D if A = Jam- 



Here H is defined by 

{(H - z)- 1 f)(x) = / dx' K 1A (z,x,x')f(x'), zeC\l,/£l 2 ( 



and this results in the following matrix-valued Schrodinger operator in L 2 
d 2 



^ _ ^ m 'dx 2 ®' 
dom(tf) = {g G L 2 (M) m | <?,</ G AC loc ( 



(2.49) 
) 2m with 

(2.50) 

(2.51) 

0"\ 
(2.52) 



s.-a. b. c. at ± 
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(-l m g" + Qg)eL 2 (R) m }, 

where "s. s.-a. b. c. at ± oo" denotes a separated self-adjoint boundary condi- 
tion at +00 and/or — oo (if any) induced by M + (z,xq) and/or M-{z,xq). Here 
Ki^i(z, x, x') is the left upper m x m submatrix of K(z,x,x'), obtained from the 
partitioning K{z,x,x') = (Kj,k(z, x, x')) 2 k=1 . Similarly, D defined by 



{{D-z)- x ^)(x) 



dx' K(z, x, x')ip(x'), z G C\R, ijieL 2 (l 
represents the following Dirac-type operator in L 2 (R) 2m , 

D=J^-- B, 

dx 

dom(D) = {(f) E L 2 (R) 2m I 4> S AC loc (R) 2m ; (J0' - Bcf)) e L 2 ( 



2 m 



(2.53) 



(2.54) 



\2m 



}■ 



where we took into account the limit point property of Dirac-type operators as 
described in Lemma 2.9, Thus, H and D are self-adjoint operators in L 2 (R) m and 
L 2 (R) 2m , respectively. 

The basic results on M(z, Xq) then read as follows. 

Theorem 2.11 (jen), ||], @, ||). Assume Hypothesis |j| and suppose 
that z e C\M, a: e R. TTien 

(i) M(z,j:o) is a matrix-valued Herglotz function of rank 2m with representation 



M(z,x ) = F(x )+ / ^(A,x )((A-z)- 1 -A(l + A 2 )- 1 ), 



whe 



F(x ) = F(x Q y 



||df7(A,x )|| C 2 m x2 m (1 + A 2 )" 1 < oo 



(2.55) 



(2.56) 



fi((A, fi], Xq) = limlim — / dv \nv(M(v + ie, xq)). 

HO elQ 7T J x+S 

(if) z £ p(T) if and only if M{z, xq) is holomorphic near z. 



(2.57) 



The fundamental uniqueness theorem for T in terms of M(z,xq) then reads as 
follows. 



Theorem 2.12 (p 



61 



116 1). Letxo S R and assume Hypothesis 



Then M{z, xq) for all z 6 C+ uniquely determines T. In particular, in the case 
of the Schrddinger operator H, M(z,xo), z G C+, uniquely determines Q(x) for 
a.e. and in the case of the Dirac-type operator D, M(z,xo), z £ C+, uniquely 

determines B(x) = ( f, 1,1 ,^, B n' 2< v\ ] for a.e. 

v ' \ -81,2(21) -Bi,i(x) J ■> 

In the context of Schrodinger operators, Theorem 2.12 represents the m atrix - 
value d extension of the celebrated Borg-Marchcnko uniqueness theorem |24|, [ 104 , 
|l05| (see also Q , [ 134 1 ) for Schrodinger operators on a half- line in the sc alar case 



m — 1. The matrix Schrodinger case has been treated by Rofe-Beketov [ 1 16 1 (as- 
suming Q to be continuous on IR) and more recently in |6(]|] (in the general case 
Q e LH$, R]) mxm for all R > 0). In the Dirac case, and for B continuous, The- 
orem 2.12 is contained in p3|. The Dirac result under the weaker Hypothesis 2.2 
is due to |pn| . In the scalar case m — 1, boundedness of B is immaterial for the 
uniqueness result to hold in Theorem 2.12 (cf. the discussion in |37]]). However, 
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for m > 2, the corresponding uniqueness result appears to have been proved only 
assuming B G Lj^ c (R) 2mx2m thus far. This is the reason why we added the local 
essentially boundedness assumption of B to Hypothesis 2.2. As soon as one suc- 
ceeds in removing this L^ c -assumption of B in Theorem 2.12| , it can be removed 
everywhere in this paper. 

The recent new results on local Borg-Marchenko uniqueness theorems for matrix- 
valued Schrodinger and Dirac-type operators on a half- line in p7], BCJ] then yield 



the following local extension of Theorem 2.12 



Theorem 2.13 ( |37|) , |60|). Suppose Bj, j = 1,2, satisfy Hypothesis |£L§ and let 
a > 0, a± > 0, and xq G R. 

(i) Denote by Mj,±(z, xq), j — 1,2 the mxm Weyl-Titchmarsh matrices associated 
with the corresponding Hamiltonian systems (2.4) on (2:0, ±oo). In the Schrodinger 
case assume that 



\\Mi t ±(z,Xo) - M 2 ,±(z,Xo)\\c" 



Q( e -2lva(z 1 ' 2 )a ± \ 



(2.58) 



along a ray with < arg(z) < tt. In the Dirac-type case suppose that 



\\M lt± (z,x ) - M 2 ,±(z,x )\\ C r, 



(9^ e -2Im(z)ai 



(2.59) 



\z\ — ^00 

along a ray with < arg(z) < ir/2 and along a ray with tt/2 < arg(z) < tt. Then, 
in either case, 



Bi(x) = B 2 {x) for a.e. x G [xq, Xq ± a±\. 



(2.60) 



(ii) Denote by Mj(z, x ), j = 1, 2 the 2mx2m Weyl-Titchmarsh matrices associated 
with the corresponding Hamiltonian systems (2.4) on K. In the Schrodinger case 
assume that 



\\Mi(z,x ) - M 2 {z,x )\\ C 2 m x2 



0(e- 2Im (* 1/2 )«) 



along a ray with < arg(z) < tt. In the Dirac-type case suppose that 
\\M 1 (z,x Q )-M 2 (z 1 x Q )\\ C 2 m ,2 m = 0(e- 2Im ^ 1/2 ) Q ) 

\z\— >oo 



(2.61) 



(2.62) 



along a ray with < arg(z) < tt/2 and along a ray with tt/2 < arg(z) < tt. Then, 
in either case, 



Bi{x) = B 2 {x) for a.e. x G [xq — a, xq + a] 



(2.63) 



Proof. In the context of Schrodinger operators, case (i) is proved in |p0| ; in the 
corresponding Dirac -type situation case (i) is proved in |37j . Case (ii) then follows 
by combining ( 2.47 ) and case(i), taking into account the asymptotic expansions 
(along any ray with < arg(z) < tt) 



M±(z,x ) = ±iz 1 / 2 I m + o(l) 

\z\ — >oo 

in the case of Schrodinger operators (cf. |56|) and 

M±(z,x ) = ±U m + o{\) 

\z\ — >oc 

in the case of Dirac-type operators (cf. |57| ) . 



(2.64) 

(2.65) 
□ 
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In the case of scalar Schrodinger operators, Theorem 2.13| (i) is due to Simon 
|134| . An alternative proof was presented in |30| , and more recently in jl6| . Most 
recently, Alexander Sakhnovich kindly informed us that his integral representation 
of the Weyl-Titchmarsh matrix in |12C[ can be used to derive asymptotic expansions 
for the Weyl-Titchmarsh matrix and its associated matrix-valued spectral function, 
and also yields a result analogous to Theorem 2.13) (i) for a certain class of canonical 
systems. Moreover, in the case of skew-adjoint Dirac-type systems, similar results 
are discussed in [ 121 1 and applied to the nonlinear Schrodinger equation on a half- 
axis. 

3. Uniqueness Results for Schrodinger and Dirac-Type Operators 

In this section we prove our uniqueness results for matrix-valued Schrodinger 
and Dirac-type operators. 

Since in the case of Schrodinger operators, B{x) is of a very special nature, 

B(x) = (~ Q ^ x) , Q(x)* = Q(x), (3.1) 

and only the mxm submatrix B\\ — —Q(x) contains information on Q(x), we will 
focu s on the m xm submatrix M\ t i{z, Xq) of M(z,xq) in ( [2.47j ) and J2.55| ), ( 2.56 ). 
By ( 2.46| ) and ( [2.51 ) one infers that the Green's matrix G(z, x, x') of H is given by 

G(z,x, x') — Ki^i{z,x,x') 

= ij) Tj \(z,x,xo)[M-(z, x ) - M + (z,aso)] - V±,i(^) x',x )*, 

x = x', z e C\R, (3.2) 



where the mxm matrices ip±^(z, x,xq) are defined in (2.44), that is, 

tp± t j(z,x,x ) = Oj{z,x,x ) + <j>j{z, x,x )M±(z, x ), j = 1,2. (3.3) 

For future purpose we abbreviate the diagonal Green's function G(z, x, x) of H by 
g(z, x), that is, 

g(z,x) = G(z,x,x), zeC\l, xeR. (3.4) 

Moreover, since "diagonal elements" of matrix- valued Herglotz functions are (lower- 
dimensional) matrix- valued Herglotz functions (cf., e.g., Q), Mi ) i(z, Xq) = g(z,x ) 
is an mxm matrix-valued Herglotz function satisfying 

Mi, i(z, x ) = g{z, x ) = Ki t i(z, x ,x ) (3.5) 

= [M^(z,x )-M + {z,xo)]- 1 (3.6) 



= F 1 , 1 (x )+ / dQ 1 , 1 (X,xo)((X-z)- 1 -Xil + X 2 )- 1 ), (3.7) 

where 

F hl (x Q )=F 1A (x Q )*, [ ||dr! 1 , 1 (A, a; o)|| Cm x m (l + A 2 )- 1 <oo (3.8) 

and 

Oi i((A, n], xq) = limlim - / dvIm(Mi i(u + ie, Xq))' (3-9) 

<510 ej.0 e J x+S 

In complete analogy to the Schrodinger case, we also introduce 

g{z,x )=M 1A (z,x ) = [M_(z,x )-M + (z,x )}-\ z e C\R (3.10) 
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in the case of Dirac-type operators D and recall that B(x) is now of the type 
'Bi,i(a;) B h2 (x 



B{X) ~ \B 2A (x) B 2 , 2 (x) 
with B ia {x)* = B 2A (x), B jtj (x)* = B jd (x), 3 = 1,2. 



(3.11) 



Next, we provide a quick derivation of the matrix Riccati-type equations satisfied 
by M±(z, x) in the Schrodinger and Dirac cases. This result is not new (see, e.g., 
|38| , f35fl , f7q] ). We only mention its proof in some detail to reach a certain level of 
completeness and since it is a fundamental ingredient in our principal uniqueness 
result, Theorem 3.7. 



Lemma 3.1. Assume Hypothesis 2.1, z G C\R, and x G K. Then in the case of 
matrix-valued Schrodinger operators H , M±(z, x) satisfies the standard Riccati-type 
equation, 

M' ± {z,x) + M±(z, x) 2 = Q{x) - zl m (3.12) 
for a.e. x € K. In the case of Dirac-type operators D, M±(z,x) satisfies 

M'±{z,x) + zM±(z,x) 2 + M±(z, x)B 2 ,2(x)M ± (z,x) + B 1 ^ 2 (x)M ± (z, x) 
+ M ± (z,x)B 2 . 1 (x) = -B hl (x)-zl m (3.13) 
for a.e. x € M. 



Proof. In order to prove ( 3.12 ), let 

'ip- t i(z,x,x ) ip+,i{z,x,x ) 
M>-,2(z,x,x ) ip + , 2 (z,x,x ) 



$>(z,x,x ) 



z G 



(3.14) 



be the fundamental system of solutions of ( 2.12 ) as defined in ( 2.44 ) and observe 
that 



ip + . 2 (z,x,x a ) = ip' +1 (z,x,x ) 



(3.15) 



in the Schrodinger operator case. Thus, any nonnormalized solutions ip±,i(z, ■ ) G 



i 2 ((±oo, c)Y 



of —ip" + Qtpi = zi/ji for z G C 
4>±,x{z,x) = ip±,i(z,x,x )C± 
for some nonsingular m x m matrices C±. Thus, 

i>±,i(z,xo) = C±, ip± A (z,Xo) = M±(z,x Q )C± 
by ( |2.15| ) and ( |2.44| ) . In particular, 

M±(z,x ) = ~ip'±^(z, x )4>± t i(z, xq)' 1 , z G C\I 



are of the type 



(3.16) 
(3.17) 
(3.18) 



is independent of the nor maliz ation chosen for ip±,i(z, xq). Varyi ng the reference 
point xq G R then yields ( 3.12 ). Similarly, in order to derive ( 3.13| ), let 



^(z,x,x ) 



ip- t i(z,x,x ) ip +A (z,x,x ) 
ip- i2 (z,x,x ) ifj +i2 (z,x,x ) 



z G C+ 



(3.19) 



be the fundamental sys tem o f solutions of ( 2.12 ) as defined in ( 2.44 ). Any other 
fundamental system of ( 2.12 ) of the form 



S&(z,x) 



$- t i(z,x) ip +A {z,x) 
ip- i2 (z,x) ib +i2 (z,x) 



z G 



(3.20) 



UNIQUENESS FOR SCHRODINGER, JACOBI, AND DIRAC-TYPE OPERATORS 15 



with ty±(z, x) = (ip±,i(z, x )i $±,2(2, x)) 1 a basis of N(z, ±00) for z G C + , necessar- 
ily must be of the form 

^>{z,x) = ^{z,x,x )G. (3.21) 

Here C is of the type 

c =( C o + 1) < 3 ' 22 > 

and C± are invertible m x m matrices. In particular, 

$-,j( z > x ) = ip-,j(z,x,%o)C-, 4> + j{z,x) = ip +J (z,x,x )C+, 3 = 1,2 (3.23) 
and thus, 

$±,i(z,Xo) = C±, i>±, 2 (z,Xo) = M±(z,x )C± (3.24) 
by ( 2.15| ) and Q2.44 ). Consequently, 

M±(z, a; ) = i>±,2(z, x )i>±,i(z, xoy 1 , zeC + (3.25) 
is independent of the normalization chosen for ip± j . Varying xq £ R one verifies 

O- □ 



Corollary 3.2. Assume Hypothesis 2.1, z £ C\R, and x £ R. Moreover, define 

g(z,x) = [M_{z 7 x)-M + {z 7 x)}- 1 , z 6 C\M, iff (3.26) 

as m (3.4) and ( |3 .10 ). Then in the case of matrix-valued Schrddinger operators, 
the diagonal Green's matrix g(z, x) of H , satisfies 

g'(z,x) = g{z,x)Mj r {z,x) + M_(z,x)g(z,x), (3.27) 

g\z,x)Tl m =9{z,x)M ± {z,x) + M ± {z,x)g{z,x) (3.28) 

for all In the case of Dirac-type operators D, g{z, x) satisfies 

g (z, x) — zg(z, x)M + (z, x) + zM-(z, x)g(z, x) — g{z, x)M±(z, x)Bx,t(x) 

- Bi t i{x)M T {z,x)g{z,x) + g(z, x)B ia {x) + B 1 ^ 2 {x)g(z,x) (3.29) 

for a.e. x S R. 



Proof. In order to prove (3.27) one only needs to observe 

(d/dx)(g(z,xy 1 ) = -g{z,xY 1 g\z,x) giz.x)" 1 

= M'_(z,x) - M' + (z,x) 

= M + (z,x) 2 - M^(z,x) 2 

= M + (z,x)[M + (z,x) - M-(z,x)] 

+ [M + (z,x) - M_(z,x)]M_(z,x) 

= -M + (z,x)g(z,xy 1 - g(z,x)~ 1 M_(z,x) (3.30) 

for a.e. x € R, using (3.12) and ( |3.26| ). By continuity of the right-hand side of 
(3.3C) with respect to x, (3.30) extends to all x £ R. Equation (3.30) then proves 
(3.27), and using (3.26) again, (3.27) proves (3.28). Similarly, to prove ( 3.29| ), one 
computes 

g'(z,x) = ~[M-(z,x) - M + {z,x)]~ 1 [M'_{z,x) - M' + (z,x)} x 
x [M-(z,x) - M+(z,x)]~ l 
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= ~[M-(z,x) - A// + (z,x)]~ 1 [zAf + (z,2;) 2 - zM^(z,x) 2 

- M + (z, x)B lA (x)M + (z, x) + M-(z, x)B x ,i(x)M-(z, x) 
+ B li2 (x)M + (z,x) - B h2 (x)M_(z,x) + M+(z,x)B 1 . 2 (x) 

- M-(z,x)Bi t2 (x)][M-(z,x) - M+(z,x)]~ l 
= zg{z, x)M + (z, x) + zM-(z, x)g(z, x) 

+ [M-(z,x) - M + (z,x)r 1 [M + (z,x)B 1 s(x)M + (z,x) 
+ M±(z,x)Bi A (x)M T (z,x) - M±(z, x)Bi A (x)M T (z, x) 
+ M_ (z, x)B-i > i(x)M- (z, x)} [M_ (z, x) - M + (z, x)}^ 1 
+ g(z, x)B li2 (x) + Bx t2 (x)g(z, x) 
— zg(z, x)M + (z, x) + zM-(z, x)g(z, x) — g(z, x)M±(z, x)Bi_i(x) 

- Bx t x(x)M T (z,x)g(z,x) + g(z, x)B 1<2 (x) + Bx >2 {x)g{z, x) (3.31) 
for a.e. 16R, using |DJ) and (|^). □ 

Denoting by spec(A) the spectrum (i.e., the set of eigenvalues) of an n x n matrix 
A G C™ x ", we next recall a well-known lemma concerning the solution of Sylvester's 
matrix equation. We also recall that A is said to be accretive if Re((x, Ax)c n ) > 
for all x G C™. Finally, we recall that 1 1 1 • 1 1 1 is said to be a unitarily invariant norm 
on C nxn if |||E/AV||| = IPHI for all U,V G C" xn unitary and all A G C" xn . 

Lemma 3.3 (c.f. [pj, Theorems VII.2.1 and VII.2.4 and [||). Assume A, B, C G 
C" x " for some n€N. 

(i) Suppose that 

spec(^) n spec(B) = 0. (3.32) 

Then the equation 

AX-XB = C (3.33) 

has a unique solution X G C nx ™. Moreover, suppose T is a closed counterclockwise 
oriented Jordan contour such that spec(^4) has winding number +1 and sj> ec(B) 
has winding number with respect to T. Then an explicit solution for X in (3.33) 
is provided by 

X = (2m)- 1 j d({A- C) _1 C(S - (3.34) 

(ii) Assume that A — (5/2)I n and —B — (8/2)I n are accretive. Then (3.33) has a 
unique solution X G C nx ™ satisfying 

\\\X\\\ <&-\\\C\\\ (3.35) 
for any unitarily invariant norm \ \\ ■ \ \\ on C nxn . 



Without going into details, we note that Lemma 3.3 (i) remains valid for bounded 
operators on a Banach space and also extends to certain cases where A and B 
are unbounded (cf., e.g., Jll5| ). Similarly, Lemma |3.3| (ii) is proven for maximal 
accretive operators A — (6/2)1, —B — (5/2)1 in a Hilbert space in p3| . 

Next, we recall one more auxiliary result to prove the principal uniqueness results 
of this section, Theorems 3.7 and |3.8| . Although we are only applying the following 
result in the special case of to x m matrices, we state the corresponding facts in 
the infinite-dimensional context since, in contrast to Lemma 3.3, this appears to be 
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a less familiar result. In the following let ri be a complex separable Hilbert space 
with scalar product (•, -)n (and norm || • \\n), an d denote by B(7i) the Banach space 
of bounded linear operators on ri. 

Lemma 3.4 (c.f. p6| , Lemma 2.3). Let M(z) be an operator-valued Herglotz func- 
tion with values in B(Tt) (i.e., M is analytic in C+ with Im(Af(z)) > 0). // 
(Im(M(z ))) _1 £ B(H) for some z G C+, then (Im(M(z)))- 1 G B(H) for all 
z G C+ and 

spec(M(z)) C C+ and hence spec(M(z)) nK = 0. (3.36) 

Proof. Suppose there is a z\ G C+ such that (Im(M(zi))) _1 ^ B(H). Then there 
is a sequence {e„}„ C H, \\e n \\n = 1, n G N, such that 

lim \\Im(M(z 1 ))e n \\ H = 0. (3.37) 

n — >oo 

Any Herglotz function M(z) with values in B(Ti) admits the representation 

M{z) =A + Bz + R 1/2 {L K + zL)(L - z)- x R l ' 2 \ H , (3.38) 
where A = A* G B(H), < B G B(ft), /C is a Hilbert space, K, C H, L — L* is a 



self-adjoint operator in /C, and < i? G B(JC) with i?|jcew = 0. By ( 3.38 ), 

Im(M(z)) - lm{z)(B + R 1 ' 2 ^ + L 2 )({L - Re(z)) 2 + (Im^)) 2 )-^ 1 / 2 ^) 

= Im(z)(B + R 1 / 2 C(z)R 1 / 2 \ n ), (3.39) 

where C(z) = (/ JC + L 2 )((L-Re(z)) 2 + (Im(z)) 2 )- 1 > is invertible in K, C(z)- 1 G 
B(/C), z G C+. Thus, ( pT37| ) together with ( pS9[ ) imply 



lirn ((en, Be„)« + (i? 1/2 1 w e„, C^R 1 ^ \ H e n ) n ) = (3.40) 
and hence lim^oo \\B 1 / 2 e n \\ H = Hm n _» 00 ||Se„|| w = lim^oo ||.R 1/2 |HeJ|« = 0. 



Applying (3.39) again one infers limn-^ ||Im(M(z))e„||-H = for all z G C+, con- 
tradicting the hypothesis (Im(-M(zo))) _1 G B(7Y). 

Since Im(M(z)), z G C+ is invertible, the numerical range W z of the operator 

M(z) 

W z = {(M(z)e, e) H \ e G W, ||e||„ = 1} (3.41) 

is a subset of the half-plane {( G C+|Im(C) > ||(Im(Af(z))- 1 )||g ( 1 H) }. By the 
spectral inclusion theorem (see, e.g., |33|, Theoreml.2-1]), the spectrum of M(z) is 



contained in the closure of the numerical range W z and hence (3.3E) holds. □ 

Moreover, we will also apply the following elementary result. 

Lemma 3.5. Let A,Be C" xn for some n G N and 

p||c»*» < 1/2, ||B|| c -xn < 1/2 (3.42) 

for some matrix norm j • | on C nxn . Then the matrix-valued Riccati equation 

Y = YAY + B (3.43) 

has a unique solution X G C nxn nK, where K = {X G C nxn | |]X|| C nxn < 1} is 
the closed unit ball in C nxn . Moreover, the solution X is given by the convergent 
series 



X = ^B k , (3.44) 



k=0 
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where 

B = B, B k+1 = B k AB k +B, ke N . (3.45) 
Proof. Introducing the map 

F : C nxn — > C nxn , F(X) = XAX + B, (3.46) 

the estimate 

||F(X)|| c „xn = \\B + XAX\\ C ,^ < ||fl|| c »x« + ||X||2 nX „||4|| c „xn < 1, XeK, 

(3.47) 



combined with ( 3.42 ) shows that F maps K into itself. Moreover, using 
\\F{X) - F(Y)\\ C nxn = \\XAX - YAY\\ c ^ n 
= \\XAX - XAY + XAY - YAY\\ c ^ n 

< (||a:||c»x. + ||y|| c -x») \\a\\ c ^\\x - y[[ c „ x „ 

<2||A|| c „x„[[X-F|| c „xn, ITel, (3.48) 

( [3.42] ) implies that F is a strict contraction. An application of Banach's fixed point 
theorem then proves the existence of a unique fixed point X £ K . Moreover, the 
fixed point X can be obtained as the limit of the sequence {i?fc}fc e N, where 

B k+1 =F(B k ), fcGN , (3.49) 

and any fixed Bq £ K, 

X = lim F(B k ). (3.50) 

k — >oc 

Taking Bq = B proves representations ( [3.44 ) and ( 3.45| ), which in turn completes 
the proof since any fixed point of F is a solution of the matrix-valued Riccati 
equation ( 3.43| ). □ 



Corollary 3.6. Assume Aj, Bj, j — 1,2, satisfy the assumptions of Lemma 3 
and Xj, j = 1,2, are the corresponding (unique) solutions of the matrix-valued 
Riccati equations 

Y = YA j Y + B j , j = l,2 (3.51) 
in the closed unit ball K . Then 

Ai - A 2 c «*» < — — — . (3.52) 

1 - 2mmj = i,2 H4j'IIc™x™ 

Proof. For Xj £ K , j = 1, 2 the estimate 

ll-Xj. — ^Hc-xn = \\Bi — B 2 + X\A\X\ — A"2A 2 A"2||c"X" 

< - S 2 || c „x„ + HXi^Xi - X 2 A 2 X 2 \\ C n*« 

= H-Bl — _B 2 || C ,ixn + ||XlAlXl — X<lA\X>l\<£ny.n 
+ ||A" 2 AiA~ 2 — -X" 2 j4 2 A~ 2 ||(>Xn 

< - B 2 || C „x„ + 2||A X [[ C nXn||X 1 - X 2 || C „x„ 

+ ||x 2 ||| nxn p 1 -A 2 || c „ x „ 

< 11^ - A 2 || c „x„ + - B 2 || c „x„ + 2||A 1 || c „x„||X 1 - X 2 || c „x„ (3.53) 
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shows that 

< 



1!) 



\\Xi — X2||oxr. 
Analogously, one proves that 

\\Xi - A 2 ||c™xr. 



< 



\A, - A 2 \\ c ^ n + WB, - B : 
l-2||i4i||c»x» 

-\\B 1 -B 2 \ 



\A X -A 



211C" 



1-2IM 



Combining (|||J) and ^5^) yields (fl||). 



(3.54) 

(3.55) 
□ 



After these preparations we are able to prove the following fundamental unique- 
ness result. 



Theorem 3.7. Assume Hypothesis 2.i and let xo G M. 

(i) In the case of a matrix-valued Schrddinger operator H, g(z,xa) and g'(z,xo) 
for all z G C+, uniquely determine the matrix-valued m x m potential Q(x) for a.e. 
x g E. 

(ii) In the case of a Dirac-type operator D assume xq is a point of Lebesgue con- 
tinuity of B and define g'(z,xo) by the right-hand side of (3.29). Then g(z,xo), 
g'{z,xo) for all z G C+, and the m x m coefficients Bi^i(xq), Bi^(xq) uniquely 

determine the 2m x 2m matrix B(x) = ( n 1 ' 1 ^ B a' 2 ^\ I for a.e. x £ IR. 

Proof. Since M(z, xq) uniquely determines T (i.e., H or D) by Theorem 2.12 , it suf- 
fices to show that g(z, xq) and g'(z, xq) uniquely determine M(z, Xq) in case (i) and 
similarly, g(z,Xo), g'(z,Xo), and Bi t i(xo), -61,2 (#0) un iquel y dete rmine M(z, xp) in 
case (ii). We start with the Schrodinger case (i). By ( 2.36| ), (|3.5[), and (^?|), 

ha(g(z, x )) > for all z G C+ (3.56) 

and hence 

spec(g(z,x )) n spec(— g(z, xq)) = for all z e 

B 



by 



(3.57) 

g(z, x ) and C = g'(z, x ) T In) 



36). Applying Lemma 3.3 (with A 
to ( p. 28 ) (with x = xq) then shows that M±(z,xo) are uniquely determined by 
](z,Xq) and g'(z,x ). 

In the corresponding Dirac case (ii) one rewrites ( |3.29 ) (with x — xq) in the form 

g(z,x )~ 1 [I m - z~ x Bi t x(x )]M + (z, x ) + M + (z,x )[I m - z~ 1 Bis(x )]g(z,x Q )~ 1 

= z~ 1 \g(z,x )~ 1 g'(z,x )g(z,x )~ 1 - Bi >2 (xo)g(z,x )~ 1 - g(z,x y x Bi t2 (xo) 

- zg(z,x )~ 2 + g(z,xay 1 B li i(x a )g(z,x a y 1 }. (3.58) 



Thus, for Im(z) > sufficiently large, we can again apply Lemma p.3| (with A = 
g(z,x )~ 1 [/ TO - z^ 1 B li i(x a )}, B = -[I m - z^ 1 Bx,i{x )]g{z, xq)" 1 , and C = the 
r.-h.s. of (3.58)) since clearly 

spec(g{z,xo)~ l [I m - z^ 1 Bx tl (x )]) nspec(-[/ m - z~ 1 B 1A (x )]g(z,x y 1 ) = 

(3.59) 

for Im(z) > sufficiently large. Hence M + (z,xq) is uniquely determined for 
Im(z) > sufficiently large. Since M + (-,Xq) is analytic in C+, one concludes that 
M + (z,xo) is uniquely determined for all z G C+. Since g(z,Xo) — [M-(z,Xo) — 
M + (z,xq))~ 1 is known by hypothesis, this also uniquely determines M~(z, xq) for 
all z G C+. 



In either case an application of Theorem 2.12 completes the proof. 



□ 
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Since g(z, xq) and g'{z, xq) are analytic in C+, it suffices of course in Theorem 3.7 
to know them for an infinite set of z G C+ with an accumulation point in C+. 

In the special scalar case m = 1, the uniqueness result stated in Theorem 3.7 
is well-known for one-dimensional Schrodinger operators on R and can be inferred, 
for instance, from the results in Q, Q, (it is also explicitly stated in 
58 1). The present matrix extension with m > 1, however, appears to be new. To 
the best of our knowledge, the corresponding Dirac case has not previously been 
studied in the literature from this uniqueness point of view. 

Next we turn to a local version of Theorem 3.7, our principal new result. 



Theorem 3.8. Suppose Bj, j = 1, 2 satisfy Hypothesis \2. e 4 and let a > and xo £ 
R. Introduce gj(z,xo) = [Mj t -(z, xq) — Mj t +(z, xq)]~~ , where Mj t ±(z,xo), j = 1,2 
denote the m x m Weyl- Titchmarsh matrices associated with the corresponding 
Hamiltonian systems ( }2.4| ) on (:tq,±oo). 

(i) In the case of matrix-valued Schrodinger operators Hj, j = 1,2 suppose that 
\\ gi (z,x )-g 2 (z,x )\\ cm , m + \\g[(z,xo) - g 2 (z,x )\\ Cm , m = 0(e- 21 ™^ a ) 

\z\ — >oo 

(3.60) 

along a ray with < arg(z) < n. Then 

Qi(x) = Q2(x) for a.e. x 6 [xq — a, xq + a}. (3.61) 

(ii) In the case of Dirac-type operators assume that xo is a point of Lebesgue con- 
tinuity of Bj and define g'j(z,Xo), j — 1,2, by the right-hand side of (3.29). More- 
over, suppose that 



||fll(z,x ) - g 2 (z,x )\\ c „ 



\\g[(z,x Q ) - g' 2 (z,x )\\ c „ 



0(e- 2lm ^ a ) 



(3.62) 



along a ray with < arg(z) < tt/2 and along a ray with 7r/2 < arg(z) < 7T. In 
addition, assume that 



(B 1 (x )) 11 = (B 2 (x Q )) lv (B 1 (x )) 12 = (B 2 (x Q )) u 



Then 



(3.63) 

(3.64) 
1, and 



Bi(x) = B 2 (x) for a.e. x £ [xq — a, Xq + a]. 

Proof. In the Schrodinger context one observes from Lemmas |3.3| and 
(for x = x ) that 

M ll± (z,x ) - M 2 .±(z,x ) 

= (2iri)- 1 f d( [( gi (z,x ) - Cy^g'^z^o) Tl m ](-9i(z,x ) - C)' 1 
Jt 1i2 

- (g 2 (z,x ) - ()- 1 lg 2 (z,x )Tl m }(-g2(z,x ) - -1 ]. (3.65) 

Here L 12 is a closed and counterclockwise oriented Jordan contour which encircles 
spec(<?i(z, Xq)) U spec (g 2 (z, x Q )) with winding number +1 and spec(— gi(z, x )) U 
spec(— g 2 (z, xq)) with winding number 0, and whose existence is guaranteed by 
Lemma 3.4 (cf. ( |3.57 )). Using the asymptotic expansion (cf. pq ]) 

gj(z, xq) = (i/2)z-^ 2 I m + oiz- 1 ), < arg(z) <ir, j = 1, 2, (3.66) 

\z\ — *oo 
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(3.6C) and a standard e/3-argument yield 
\\Mi t ±(z,x ) - M 2 ,±(z,x )\\ c 



= 0(e- 2Im ( zl/2 ) 



(3.67) 



along a ray with arg(z) = 7r — e. Hence one obtains fl3.6ip b y Theorem |2T| (i) . 
In the case of Dirac-type operators one infers from (3.58) that 

5j(^2Co) -1 [^m - ^ 1 (B j (xo)) lil ]M j;+ (z,Xo) 

+ Mj >+ (z,xo)[I m - z^ 1 (Bj(x )) 1A ]gj{z,x y 1 

= z~ 1 [g j (z,x o y 1 g' j (z,xo)gj(z,x y 1 - (B j (x Q )) 1 ^(z^xo)' 1 

- g j (z,Xo)^ 1 (B j (x )) h2 - z gj (z,x y 2 + g j (z 7 x y 1 (B :j (xo)) 11 gj{z,xo)~ 1 7 

J = 1,2 (3.68) 

and then proceeds in exactly the same manner as in (^65|)-(|3~67|) replacing ( ^66|) 
by (cf. |3) 

g 3 {z, x ) = (i/2)I m + o(l), < arg(z) < tt, j = 1, 2. (3.69) 

I z\ — >OC 

□ 



For an alternative proof of (3.67) one can multiply (3.2S) (for x = xq) and (3.68) 
by (— i) and directly apply Lemma 3.3(h) taking into account ( 3.66| ) and (3.69), 
respectively. 

As mentioned in fl6(j| |, the ray with < arg(z) < tt in Theorem 3.S (i) is of no 
particular significance and can be replaced by any non-selhntersecting Jordan arc 
that tends to infinity in the sector e < arg(z) < tt — e for some < e < ir/2. The 
analogous comment applies to Theorem 3.8 (ii) 



As discussed in Marchenko's monograph [[106|, p. 193-19 5], it is possible to reduce 
certain Dirac operators D with potentials B(x) of the type ( 2.10| ) to 2mx2m matrix- 
valued Schrodinger operators by essentially taking the square of the Dirac operator 
D. This permits one to relate the M- functions of the 2m x 2m Schrodinger and 
Dirac-type operators and we next explore this connection in some detail. In order 
to indicate that B is temporarily of the special form (2.1C), we denote it by B in 
the following and assume in addition to Hypothesis 2.2 that 



B 



Bl^i #1,2 

B\ o — B\ 1 



2m X 2m 



(3.70) 



(this can be improved but is of little relevance in our context). Denoting the 
resulting Dirac operator by D, 

d 



D 



J-j B, 

ax 



one computes, using J 2 = —hm and JB(x) + B(x)J = 0, 



'2m 



dx 2 



B 2 - JB' 



(3.71) 



(3.72) 



and hence obtains the 2m x 2m matrix-valued Schrodinger operator 
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where 

Q{x) = B{xf - JB'{x) (3.74) 

B 1 , 1 (x) 2 +B 1 , 2 (.x) 2 +B' 1 2 {x) B 1 , 1 (x)B 1 , 2 (x)-Bi, 2 (x)Bi,i(x)-B' 11 (x) 



_ B la (x)B 1A (x)-B 1A (x)B li2 (x)-B' 1A (x) B 1A (x) 2 + B 1 . 2 (x) 2 -B' 12 (x) 

Since 

{H-z 2 } = (D-z)(D + z), (3.75) 

linearly independent solutions <£> + (z 2 , x, Xq), $_ (z 2 , x, Xq) of (H — z 2 )$(z 2 , x, Xq) = 
are of the type 

^ ,2 \ f tp±,l(z,X,X ) 1p±,l(-Z,X,X )\ _ to /o 



where, in accordance with (2.44), the m columns of ^f±(z, x, Xq) 



, 2 e C\R, x G 



^±,2(2, x, x ) y 

l(z,X,X ) <pl(z,X,X )\ ( In 



h{z,x,x ) 4> 2 {z,x,x ) J \Mg (z,x ) 

form a basis for the nullspace of (D — z) Here, in obvious notation, M±(z,x ) 
abbreviates the corresponding m x m half- line Weyl-Titchmarsh matrices associated 
with D, given by 

M^z.xo) = t/j±,2(z,x Q ,xo)ip± : i{z,X( h x Q y 1 = i>±,2(z,x ,x ), 

Mg(-z,x ) = M£(-z,x )*, zeC+, (3.78) 

using the normalization ip±j(z, xq,xq) — I m , j = 1,2. The corresponding 2m x 
2m half-line Weyl-Titchmarsh matrices M± (z 2 ,xo), associated with H, are then 
defined by 

M^{z 2 1 xo) = <P±(z 2 1 x 0l xo)^±(z 2 ,xo,x y\ z 2 eC+. (3.79) 
Next we explicitly compute M±(z 2 ,Xo) in terms of M±(z,Xq), M±(—z,xq). 
Lemma 3.9. Assume Hypothesis \2.<\ , (|3.70 ), and arg(z) 6 (0, 7r/2). Then 



Mi(z 2 ,x ) = (Mf (z^ok*)^^, (3.80) 
Mf (z 2 , 10)1,1 =5i, 2 (s ) 

+ 2«Af|(z,i )[M|(-2,a:o) - m|(z, aro)]"^! (-«, ap ), 
Mf (z 2 , 10)1,2 = -Bx.ifco) 

- Z [Af|(-«,a:o) + Af|(z,a;o)][M|(-z,a:o) - m|(z, xq)]" 1 , 
M| (z 2 ,aro) 2 ,i = -B hl (x ) 

-z{M°(-z,x ) + M°(z,xo)}- 1 {M°(-z,x Q )~M°(z,x )}, 
m| (z 2 ,x ) 2 . 2 = -B h2 {x ) + 2z[M£{-z, x ) - M£{z, x )]-\ 
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Proof. Combining ( [3.7£ ), 



$±(z 2 ,x ,x )= [ 5 Im b j m ), (3.81) 

\M±(z,x Q ) M±{-z,x )J 

®±{z 2 ,x ,x Q )- x (3.82) 

= ( [M£(-z )Xo ) - Mg(z,x )]Mj>(-z,x ) -[M£(-z,x ) - m£(z, xv)}- 1 ^ 
\-[M£(-z,Xo)-M£(z,x )}M£(z,x ) [M^(-z,x )-M^(z,x )]- 1 J ' 

J*' ± (z,X ,X ) = (~ Bl,l(x ) + zI m + B 1 , 2 (x )M^Z,X ) \ ^ 

yields ( |3.8C| ) after straightforward matrix computations. □ 

4. Matrix- Valued Jacobi Operators 

In this section we prove analogous uniqueness theorems for matrix-valued Jacobi 
(i.e., second-order finite difference) operators. We could consider analogous Dirac- 
type difference operators and hence introduce discrete Hamiltonian systems as in 
Section ^, but for simplicity we will confine ourselves to the Jacobi case. 

As the basic hypothesis of this section we adopt the following. 

Hypothesis 4.1. Let to G N and consider the sequences of self-adjoint to x to 
matrices 

A = {A(k)} keZ , A{k) = A{k)*, A{k) > 0, k 6 Z, 

B = {B{k)} k ei, B(k) = B(k)* , k G Z. (4.1) 

Moreover, we assume that A(k) and B{k) are uniformly bounded with respect to k, 
that is, there exists a C > 0, such that 

\\A(k)\\ Cm , m + \\B(k)\\ Cmxm <C, keZ. (4.2) 

Next, denote by S ± the shift operators in t°°{Z) m and £°°(Z) mxm , that is, 

(S ± g)(k) = g ± (k) = g(k ± 1) for g G £°°(Z) m or g G £°°{Z) mxm , k G Z. (4.3) 



Given Hypothesis |4.1| , the matrix- valued Jacobi operator H in £ 2 {Z) m is then de- 
fined by 

H = AS+ + A~S- + B, dom(H) = £ 2 (Z) m . (4.4) 



Because of hypothesis (4.2), H is a bounded symmetric operator and hence self- 



adjoint. In particular, the difference expression AS + + A~S~ + B induced by (4.4) 
is in the limit point case at ±oo. We chose to adopt ( |4.2D for simplicity only. 
Our formalism extends to unbounded Jacobi operators and to cases where the 



difference expression associated with (4.4) is in the limit circle case at +oo and/or 



— oo. We just note in passing that without assuming (4.2), the difference expression 
AS + + A~S~ + B is in the limit point case at ±co if ~J2 k ^° ka \\A(k)\\^ xm = oo (see, 



e.g., )19| Theorem VII.2.9]). 

The Green's matrix associated with H will be denoted by G(z,k,£) in the fol- 
lowing, 

G(z,k,£) = {H - z)~ l (k,£), z G C\R, k,£ e Z. (4.5) 
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Next, fix a site fco G Z and in analogy to ( 2.13| )-( pTl5 ) define mxm matrix- valued 
solutions (j)(z, fc, ko) and 9{z, fc, ko) of the equation 

A(k)ip(z, k+1) +A(k - l)ip(z, k - 1) + (B(k) - zl m )ip(z, k) = 0, z e C, k e Z, 

(4.6) 

satisfying the initial conditions 

9(z,k ,k ) = cj)(z,ko + l,k ) = I m , 4>(z,k ,k ) = 6(z,k + l,k ) = 0. (4.7) 

As in ( ^44| ), ($U$) one then introduces mxm matrix-valued Weyl solutions 
ip±(z, k, ko) associated with H defined by 

ip±(z,k,k ) = 6(z,k,k ) - <f)(z,k,ko)A(k y 1 M ± {z,k ), zeC\l, k G Z, (4.8) 

with the properties 

A(k)ip±(z,k + l,ko) + A{k - l)tp ± (z,k- l,fc ) + (B(k) - z)tp ± (z,k,k ) = 0, 

z G C\R, k G Z, (4.9) 



e C\R, 



(4.10) 



where M±(z, fco) denote the Weyl-Titchmarsh matrices associated with H. Since 
by assumption AS + + A~ S~ + B is in the limit point case at ±oo, M±(z, ko) in 
( [Q| ) are uniquely determined by the requirement ( 4.10 ). We also note that by a 
standard argument, 



det((/)(z, k, ko)) ^ for all k G Z\{fco} and z G 



(4.11) 



since otherwise one could construct a Dirichlet-type eigenvalue z G C\R for H 
restricted to the finite segment fco + lj • • • ■, k — 1 for k > ko + 1 of Z (and similarly 
for fc < fco ^ !)• Thus, introducing 

M N (z,k ) = -0(z,iV,fco) _1 6'(z,iV,fco), z G C\R, JV G Z\{fc }, (4.12) 

one can then compute M±(z, fco) by the limiting relation 

M±(z,fc )= lim M N (z,k ), z G C\R, (4.13) 

the limit being un ique since AS' + + yl^S 1- + _B is in the limit point case at ±oo. 
Alternatively, (4.8) yields 



M±(z,fc ) = -A(ko)ip±(z,k + l,fc ), z G C\I 
More generally, recalling 

det(ip±(z, fc, fc )) ^ for all k G Z and z G C\R 



by an argument analogous to that following (4.11), we now introduce 



M±(z,k) = -A(k)ip±(z,k + l,ko)ip±(z,k,ko)~ 



z G 



fc G Z. 



(4.14) 
(4.15) 
(4.16) 



One easily verifies that M±(z,fc) represent the Weyl-Titchmarsh M-matrices as- 
sociated with the reference point fc G Z. Moreover, one obtains the Riccati-type 
equation for M±(z, fc), 



M±(z,k) + A(k - l)M±(z,fc- l) _1 A(fc-l) - (B(k) - z) = 0, z G 



fc G Z 
(4.17) 



as a result of (4.9). 
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Next, we define the self-adjoint half-line Jacobi operators H±^ in £ 2 ((k 0l ±00) n 
Z) m by 

H ±M) = P ±,k HP ±,k \ iH{kot±oo)nZ)m , k G Z, (4.18) 

where P±,fc are the orthogonal projections onto the subspaces i 2 ([k , ±00) Pi Z) m . 
By inspection, H± t k satisfy Dirichlet boundary conditions at fco =p 1, 

(H + , k J)(k ) = A(k )f+(k ) + B(k )f(k Q ), 

(H_, k J)(k ) = A-(k Q )f-(k ) + B(k )f(k ), 

(H±,k f)(k) = A(k)f+(k) + A-(k)f-(k) + B(k)f(k), fe|fc ±l, (4.19) 

/ g dom(JJ±, feo ) = f((k , ±00) n Z) m 

(i.e., formally, f(k =F 1) = 0). We also introduce Weyl m-matrices m±(z,ko) 
associated with -H±,fc by 

m±(z,k ) = Qk {H±M ~ z)~ l Qk , (4.20) 
= G±, fco (z,fc ,fc ), zeC\R, fc GZ. (4.21) 

Here Qfe are orthogonal projections onto the m-dimensional subspaces ^ 2 ({fco})"\ 
fco G Z and 

G±, fco (z,M) = (fl±,*o -z)~Hk,£), zeC\E, MeZ (4.22) 

represents the Green's matrix of H±^ - 

In order to find the connection between m±(z, ko) and M±(z,ko) we briefly 
discuss the Green's matrices G± t k {z 1 k,£) and G(z, k,t) associated with -H±,fc and 
H next. 

First we recall the definition of the Wronskian W(f,g) (k) of two sequences of 
matrices /(•),. g (•) G C mxm given by 

W(f,g)(k) = f(k)A(k)g(k + l)-f(k + l)A(k)g(k), k G Z. (4.23) 



We note that for any two matrix-valued solutions (p(z, •) and ip(z, •) of (4J3) the 
Wronskian W((f(z, -)*,ip(z, -))(k) is independent of k G Z. 

In complete analogy to the scalar Jacobi case (i.e., m = 1) one verifies, 



G _, fe0 (z, mh ; T , . , 1 : r r : . : (4-25) 



r , .-^(^/c^o-l^fco-l)- 1 ^,^^-!)*, ^<fc, ( , 0A , 

(2, fc, fc + l)A(k + 1)- fe + 1)*, * < fc, 
^-(*, fc, fc + l)A(k Q + l)^(j)(z, I, fc + 1)*, I > fc, 

z g C\R, kJeZ. 

Similarly, using the fact that 

ip+(z, fc, fc )[M_(z, fc ) - Al+(z, fco)]~V-(2, fc, fco)* 
= -0_(z,fc,fc o )[Af_(z,fc o ) - M + (z,k Q )]^ 1 ip + (z,k 1 k ) 
= [M_(z,fc) -M+(z,fc)]" 1 , fcGZ, (4.26) 

and 

i4(A)^ + («, fc + 1, fco)[M_(z, fc ) - M+(z, fc )p V-(*> fc, fco)* (4-27) 
- A(k)i^-(z,k + l,k )[M-{z,k ) -M + (z,*b)]~V+(^*.*o)* = A™ fc e Z, 
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one verifies 



G{z,k,t) 



iP+{z, k, k )[M_(z, fc ) - M+{z, ko)]- 1 ^-^, £, k y, i < k, 



(4.28) 



% i/>- (z, k, k ) [M- (z,k )-M + {z, k )] - (z, i,k )*, £>k, 

z e C\R, kjGZ. 

Using ( f4~9[ ), ( |4.17| ), and (|4.2l| ), one inf ers t hat the Weyl-Titchmarsh ma- 



trices M±(z,ko) intr oduc ed by ( 4.1 3| ) (resp. ( 4.14 )) and the Weyl m-functions 
m±(z, ko) defined in ( 4.20 ) are related by 

M+{z,k ) = -m + (z,/c ) _1 -zI m + B(k ) (4.29) 

and 

M-(z,k ) = m_(z,fco) _1 . (4.30) 

As in the continuous case, ±M±(z, ko) are Herglotz matrices satisfying 

±lm(M±(z,k )) >0, zeC+, M±(z,k ) = M±(z,k )*. (4.31) 

Next, we introduce the 2m x 2m Weyl- Titchm arsh matrix M(z, ko) associated 
with the Jacobi operator H in l 2 (T) m as in ( 2.47 ) (with xo G M replaced by ko G Z, 
etc.) by 

M{z,k ) = (M^,{z,k )) jjl=12 , z G C\R, k G Z, 

where 

M hl {z, k Q ) = [Af_(z, fc ) - M+(z, k Q )}-\ 

M h2 {z,k ) = 2- 1 [M_(z,ko) - M+frkoT^M-frko) + M + (z,k )], 
M 2A {z,k ) = 2- 1 [M_(z,ko) + M + (z 1 ko)][M-(z,k Q ) - M+(z, ko)]- 1 , 
M 2 . 2 (z, k ) = M±(z, k )[M-(z, k Q ) - M + {z, k Q )]- l M T {z, k ). 



(4.32) 

(4.33) 
(4.34) 
(4.35) 
(4.36) 



Remark 4.2. We note that the Weyl-Titchmarsh matrix M(z, ko) is related to the 
Green's matrix associated with the Jacobi operator H by 









-A(ko\ 



M(z,k ) 









-A(k ] 





1 1 1 J 







M(z,k ) - 
where 

M{z,k ) 

With M(z, k ) defined in Q4.32ft , the analog of Theorem |2.12| holds 



G{z,k 0l ko) G{z,k ,k + l) 
G(z,k + l,k ) G(z,k + l,k + l) 



(4.37) 



(4.38) 



Theorem 4.3. Assume Hypothesis 4-1 om,d let k G Z. Then the 2m x 2m Weyl- 
Titchmarsh matrix M{z, ko) for all z G C + uniquely determines the Jacobi operator 
H and hence A — {A(k)} kez and B = {B(k)}ki£z- 



Perhaps the simplest way to prove Theorem 4.3 is to reduce it to knowledge of 
M±(z, ko), and hence by ( 4.29] ) and (4.30) to that of m±(z, ko) for all z G C+. Using 
the standard construction of orthogonal matrix- valued polynomials with respect to 
the normalized measure in the Herglotz representation of m±(z, ko), 



m±{z,ko) = / dil±(X, k a ) (A - z) 1 , z G C+, 



dQ±(X, ko 



(4.39) 
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allows one to reconstruct A{k), B(k), k G [fco, ±oo)nZ from the measures dfl±(X, fco) 
(cf., e.g., |1 Section VII.2.8]). More precisely, 



A(k ±k)= / \p ±tk (x,k )dn ± (x,ko)p±, k+ i(x,k y, 

JR 

B(k ±k)= [ XP ±>k (X,k )dn ± (X,k )P ±>k (X,k )*, keN , (4.40) 

JR 

where {P±,k{X, fco)}fceN is an orthonormal system of matrix-valued polynomials 
with respect to the spectral measure d£l±(X, fco), with P± t k{z, fco) of degree fc in z, 
P±,o{z,ko) — I m . One verifies, 

P+, fe (z,fco) = </>(z,ko + k,ko-l), (4.41) 
P_ jfc (z,fc ) = 9(z,k - fc,fc ), fc £ {-1}UN , fc E Z, z e C, (4.42) 



with 4>(z,k,ko) and 9(z,k,ko) defined in (16), (4.7). 

In the following we denote by [x\ £ Z the largest integer less or equal to x € M. 
The local analog of Theorem [D^ then reads as follows. 

Theorem 4.4. Suppose < e < it. Let Hj, j = 1,2 be two Jacobi operators 
in £ 2 (Z) m , with coefficients {Aj t k}ktz, {Bj t k}kei satisfying Hypothesis J^.l, and 
denote by Mj j ±(z,ko) the mx m Weyl-Titchmarsh matrices associated with Hj, 
j = 1,2. If there exists an N eN. N > 2, such that 

||Ml,_(j S) * )-M 2l _(z,to)|| C m*m , = 0(\ Z \ 



-N\ 



along a ray with arg(z) = tt — e, then 

Ai(k - fc) = A 2 (k - fc), fc = 

Bi(Ao-fc) = S 2 (*d-A), * = 
If there exists an N G N. N > 2, such that 

||M lj+ (z, fc ) - M 2 . + (z, fc )|| c ™- 



0, ... , W2J - 1, 
0, ...,L(iV-l)/2j. 



0(\z\- N ) 



along a ray with arg(z) = tt — e and 

B 1 {k ) = B 2 {k ), 

then 

Ai(*o + fc) = A 2 (k + fc), fc - 0, ... , LA72J - 1, 
B 1 (k + k) = B 2 (k + k), fc = 0, ... , L(AT - 1)/2J. 

Proof. Using (4.29), ( 4. 30] ), and the asymptotic representations 



m jt ±(z,k ) 

one infers 

mi ±(z, fc ) - m 2 .±(z, fc ) 



0(\z\ 



3 = 1,2, 



(4.43) 

(4.44) 
(4.45) 

(4.46) 

(4.47) 
(4.48) 



|z|^-oo 



O(\z\~ 2 (M h± (z,k ) - M 2 . ± (z,k ))). (4.49) 



Here assumption (4.46) is explicitly needed in the case of the +-sign in ( 4.49| ). By 
hypothesis fl4.43|), ^4~49|) implies 



m h ±(z,k )-m 2 ,±{z,k ) = 0(\z\-( N+ V), 

\z\ — >-oo 



(4.50) 
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and hence 



dn h± (X,k )X r 



dfl 2 . ± {X,k )X n 1 n = 0,l, 



,N. 



(4.51) 



Introducing the orthogonal polynomials \Pj.±-k(X, fco)} fcgN n associated with the 
normalized spectral measures dflj t ± (A, fc ), j = 1,2, ( 4.5l| ) implies that 



/ APi, ±;fe (A, fco)dn ll± (A, fc ) Pi,± :fc +i (A, fc )* 

JR 

AP2,± ; fc(A, fco)d0 2 ,±(A, fc ) P 2 ,±-k+i(X, fc )* 



0, LAT/2J - 1, (4.52) 



and 



J XPi ±. k (X, fc )dO 1)± (A, fc ) A, ±;fc (A, k )* 
AP 2 , ±; fe(A, fco)dn 2 ,±(A, fco) P 2 , ±; fe(A, k y 



which proves (4.44) using (4.40) 



fc = 0,..., L(iV - 1)/2J, (4.53) 



□ 



Of course, one can permit e = and e = 7r in Theorem |4.4| (and in Theorem 4.S 
below) as long as \z\ > max J= i j2 \\Hj\\- 

Remark 4.5. The apparent asymmetry with respect to the additional condition 
( [4.46| ) in connection with Mj, + (z, fco) is due to the fact that for fixed fco £ Z, the 
map 



(m + (z,k ),B(k )) h-> M + (z,ko), ze 



(4.54) 



given by (4.29) is not injective and hence knowledge of M + (z, fco) does not determine 
m + (z,fco) uniquely This follows from the explicit cancellation of -B(fco) in (4.29) 
as \z\ — > oo in sharp contrast to the asymptotic behavior of M_(z, fco) in (4.30), as 
described explicitly in ( 4.59| ) and (4.6C) below. 



Given these preliminaries and introducing the diagonal Green's matrix by 

g(z,k) = G(z,k,k), zeC\l, fceZ, (4.55) 



we can now formulate the analog of the uniqueness result, Theorem 3.7 , for Jacobi 
operators. 



Theorem 4.6. Assume Hypothesis and let fco G Z. Then any of the following 
three sets of data 

(i) g(z, fco) and G{z, fco, fco + 1) for all z G C+; 

(ii) g{z, fco) and [G(z, fco, fco + 1) + G(z, fco + 1, fco)] for all z G C+; 
(Hi) g(z, fco), g(z, fco + 1) for all z G C+ and A(ko); 

uniquely determines the matrix-valued Jacobi operator H and hence A = {A(k)}k<£Z 
and B = {B{k)} k£Z . 



Proof. Case (i). A direct computation utilizing formulas ( 4.28 ), ( [4.55] ) shows that 
g(z, fco) = [M-(z,ko) - M + (z,k )]-\ (4.56) 
G(z, fc , fco + 1) = -[AM*, fc ) - M+{z, fc )]^ 1 M + (z, ko)A(ko)- 1 . (4.57) 
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In addition, the matrix A(ko) can be read off from the asymptotic expansion of 
G{z, fco, fco + I, z) as z — > oo by a simple Neumann series-type argument. Indeed, 
by d420| ) 

m±(z, A: ) - -z- 1 /^ - z- 2 B(k ) - z- 3 (A(k ) 2 + B(k ) 2 ) + 0{\z\- A ) (4.58) 
and hence by ([0|) and ( |4.30| ) 

M+{z,k ) = -z^Aiko) 2 +Q(\z\- 2 ), (4.59) 



M_(z,k ) 



-zl m + B(k ) + z^Aiko) 2 + 0{\z\- 2 ) 



(4.60) 



along any r ay wit h arg(z) = 7r — e, < £ < 7r. Therefore, combining ( 4.56| ), ( 4.57 ), 
(||§), and (l4~6(il ) one infers that 



G(z,fc ,fc + 1) = -z- 2 A(k a )+0(\z\- 3 ), k e 

\z\ — >oo 



(4.61) 



along any ray with arg(z) = it — e, < e < ir, in complete analogy with the 
corre sponding scalar case ( see, e .g., [135], [136, Sect. 6.1]). Given A(ko) and using 
( i.5t ) one can easily solve ( 4.57 ) for M±(z,ho), 

M+(z, fc ) = -g(z, ko^Giz, fc , fc + l)A(fco), (4.62) 

M_(z,k ) = M + (z,k )+g(z,k Q )- 1 

= giz^o)- 1 ^ - G(z,fc ,fc + l)A(fco)), (4.63) 

which in turn determines A(k), B{k) for all fc 6 Z by Theorem fi~g| . 
Case (ii). The data (ii) are perhaps the closest discrete analog of the continuous 
version in Theorem 3.7. As before (cf. (4.61)), one can determine A(ko) from the 
asymptotic expansion of G(z, fco, fco + 1) + G{z, fco + 1, fco), 

G(z,k ,k + l) + G(z,k + l,ko) = ~2z- 2 A(k ) + O{\z\- 3 ), fc g Z. (4.64) 



Next, using (4.28), one observes that M + (z,fco) is the unique solution of the fol- 
lowing Sylvester-type matrix equation 

G(z, fc + 1, fc ) + G{z, fc , fco + 1) 

= -A(fc )- 1 A/+(z, k )g(z, fc ) - g(z, k )M+(z, fco)A(fco)- 1 (4.65) 

(cf. (3.33)) by Lemmas |3.3| and |3.4[ D etermining Af_(z, fco) from ( 4.63] ) completes 

the proof of case (ii) by Theorem 4.3 as in case (i). 

Case (hi). Since A(ko) is known, the data (iii), ( [4.28| ), and (4.55) imply 

g(z, fc + 1) = A(fc )- 1 A/+(z, fc )[M_(z, fc ) - M+(z, fc )]~ 1 Af_(z, fco)A(fco)- 1 . 

(4.66) 



Together with formula (4.56) for g(z,ko) this leads to the following equation for 

M+(z,k ), 

M+{z, fc ) + M+(z, ko)g{z, k )M + (z, fc ) = A{k )g{z, fc + l)A(k ). (4.67) 

(4.68) 



Using the asymptotics ( 4.59 ), ( 4.6C ) and the representation 
g(z,k) = -z- 1 I m + 0{\z\- 2 ), fceZ, 
one infers 

||M + (z,fc )|| Cm x m < 1, \\A(k )g(z,k + l)A(k Q )\\ c „ 



<2' 
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and \\g(z,k ) 



1 

<2 



(4.69) 



for Im(z) > sufficiently large. Thus, by Lemma 3.5, M+(z, ko) is uniquely de- 
termined by the data (iii) for lm(z) > large enough and hence for all z G C+ 
by analytic continuation. The Weyl matrix M-(z, ko) is then determined by (cf. 

(B) 

M_{z,k ) = M + {z,k )+g(z,k )-\ (4.70) 
completing the proof. □ 

Remark 4.7. We n ote that M + (z, ko) can explicitly be represented by the contour 
integral (cf. (HH)) 

M+(z,k ) 



(2ni) 



1 1 d((g(z,k Q ) - CAiko)- 1 )- 1 ^,^^ + 1) + G(z,k + l,ko)] x 
T 

\ — 1\ — l 



x (g(z,ko) + CA(k )- 1 )- 1 . (4.71) 

Here T is a closed clockwise oriented Jordan contour such that spec(A(k )g(z 7 k )) 
has winding number +1 and spec(— g{z, ko)A(k Q )) has winding number with re- 
spect to F. Without loss of generality we may assume that T C C+. The existence 
of such a contour T follows from the observation that 

spec(<?(z, k a )A(ko)) = spec(A(fc ) 1/2 g{z, k )A(k ) 1/2 ) C C+, z G C+ (4.72) 

by Lemma 



The special scalar case m = 1 of Theorem 4.6 is known and has recently been 



derived in |135j (see also |136| , Sect. 2.7]. Condition (iii) in Theorem 4.6 is specific 
to the Jacobi case. In the corresponding Schrodinger case, the corresponding set of 
data does not even uniquely determine the potential in the scalar case m = 1 (see, 
e.g., |57|, pi). 



The local analog of Theorem 3.S in the discrete context then reads as follows. 

Theorem 4.8. Suppose < e < n. Let Hj, j = 1,2 be two Jacobi operators 
in £ 2 (Z) m , with coefficients Aj — {Aj,k}kez> Bj = {Bj ^kez satisfying Hypoth- 
esis 4-1, arid denote by Gj{z,k,£), gj{z,k) — Gj(z,k,k) (k,£ G Z) the Green's 
matrices associated with Hj, j = 1,2. Moreover, assume there exists an N G N, 
N > 2, such that one of the following conditions holds. 

(0 

\\gi(z,k ) - g 2 (z,k )\\c™*™ + \\G 1 (z,k ,k + 1) - G 2 (z,k Q ,k + l)\\ c ™xm 



= 0{\z\-^); 

\z\— >oc 



(4.73) 



iii) 



\\gi(z,k ) - g 2 {z, feo)||c""« m + \\Gi{z, fc , k + 1) - G 2 (z,k , k + 1) 
+ Gi(z, k + 1, ko) - G 2 (z, k + 1, k )\\ cmxm = OQz\-W +s >); 

\z\ — >oo 

{iii) A\{ko) = A 2 (ko) is given and 

\\gi(z,k ) - g 2 (z,k )\\c™x™ + hl(z, k + 1) - g 2 (z, fc + l)|| C mxm = 

1*1- 



(4.74) 

0(|z|~W) 
(4.75) 
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along a ray with arg(z) = 7r — s. Then 

Ax(ko ±k) = A 2 (k ±k), k= 0, ... , LA72J - 1, 

B 1 (k ±k) = B 2 (k ±k), k = 0, ...,L(iV-l)/2j. (4.76) 

Proof. Case (i). Since N > 2 in fl4.73j), (J4.61[ ) implies that Ai(k ) = A 2 (k ) and 
we denote the latter by A(ko). By ( (4.62 ) the m x m Weyl-Titchmarsh matrices 
Mj :+ (z, ko) associated with Hj, j = 1, 2 can be represented as 

M h+ (z,k ) = ~g J {z,k )- 1 G J (z,k ,k + l)A{k ), j = 1,2 (4.77) 

and they admit the estimate 

||Mi, + (z, fc ) - M 2 , + (2, k )\\ C mxm 

< \\A(ko)\\ c ™x™\\gi(z, fco) _1 Gi(z, fc , h + 1) - .92(z,/co) _1 G 2 (2:,A:o,fc + l)\\ C mx m 

< ||A(fc )||c'"x m ||(7i(z,fco)- 1 || Cm x m ||G'i(z,fco,A:o + 1) - G 2 (z, k , k + l)|| c — 
+ P(fc )|| c ™x m || 3 i(z,fc )- 1 -sa^Ao) -1 !!^"' ||Ga(«,Ad,Ao + l)||c»x» 

< ||A(fc )||c—ll3i(^fco)" 1 ||c—l|Gi(z,fco,fco + l) -62(2,^^0 + l)|| Cmxm 
+ ||A(/!o)||cmxm ||5i(z, fc )" 1 || cm x m || 5 i(z,/j ) -fir 2 (2,feo)||c"»x"» x 

x || 52 (z,A : o)- 1 || cm x m ||G 2 (z,A : o,fco + l)|| C ™x m . (4.78) 
Using the asymptotics 

= -z^Ira + Odzl- 2 ), J = 1,2, (4.79) 

I z I — >oo 

Gj(z, feo,A:o + l) , = -z~ 2 A(k ) + O(\z\- 3 ), j = 1,2, (4.80) 

\z\ — *>oo 

the estimate ( 4.7S| ) combined with ( 4.73| ) proves 

||M 1)+ (*,Ao)-M 2l+ (z J *o)||c»x» i = OQzl-W+V). 



Since 



M jt -(z, k ) = Mj,+ (z, k Q ) + gj (z, ko) \ 3 = 1,2, 
and CpLSlb yield 

||M li _(z,fc )-M 2 ,_(«,A )|| C mKm = 0(\z\- N ). 



(4.44) and (4.83) yield 



(4.81) 

(4.82) 
(4.83) 

(4.84) 



S 1 (fe )=S 2 (fe ). 

The asymptotic expansions (4.81) and (4.83) together with Q4.84) ) then prove (4.76) 
by Theorem 4.4 . 

Case (ii). Since TV > 2 in (|7rJ), ( |l6i| ) again implies that A\(ko) = A 2 (ko) and we 
denote the latter by A(ko). Given e £ (0,7r), let 7 be a closed clockwise oriented 
Jordan contour 



7 C {z e C | Re(z) > 0} 
such that spec(v4(fco)) has winding number +1 and 

e 4E 7 C C+. 



(4.85) 
(4.86) 
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Next, given z £ C+, arg(z) = 7r — e, introduce the contour T(z) in the upper 
half-plane by 

T(z) = -z- 1 ^. (4.87) 



Using ( 4.7£ ) one concludes that spec(A(k ) 1 ^ 2 g 3 ;(z, ko)A(ko) 1 / 2 ) has winding num- 
ber + 1 with respect to the contour T(z) if z lies on the ray arg(z) = tt — e for 
\z\ sufficiently large. For such z the matrix-valued functions Mj t +(z, ko) admit the 
representation (cf. (4.71)) 

Mj )+ («,fco) = (2m)- 1 <t dCG i (2 ! fc ,C)[Gj(^,A;o ! A;o + l) + Gj(^fco + l,feo)] x 



r( z ) 



Gj(z,fe ,-C), J = 1,2, 



where 



Gj(z,ko,0 =^(fco) 1/2 (^(fco) 1/2 5 J (^fco)v4(fc ) 1/2 -C)- 1 ^(fco) 1/2 , 

j = 1,2, cer(z). 

Introducing 

fj(z,k ) = 9j(z, k ) + z _1 im, i = U 2 , 
one infers by (4.79) that 



/,-(*, fco) , = 0(|z|- 2 ), j = 1,2. 



(4.88) 
(4.89) 

(4.90) 
(4.91) 



The asymptotics 



sup Wi-z-'Aiko) + A(fc ) 1/2 /,(z, fco)A(fc ) 1/2 T O^llc-x- 
Cer( z ) 

= |z| sup ||(^(fco)-^(fc ) 1/2 / 3 (^A : o)A(fc ) 1 / 2 ±z0^ 1 ||c™ 
Cer( z ) 

= |z| sup ||(A(fc ) T _1 ||c-x m x 

CG7 

sup||(/ m -z(A(fc )TO" 1 ^(fco) 1/2 /,(^A ; o)A(fco) 1/2 )- 1 || c , 



C£7 



0(1*1) 



(4.92) 



|z| — >oc 

along the ray arg(z) = 7r — e, and the representation 

Gj(z, k , Q = A(fc ) 1/2 (-^ 1 A(fc ) + A{ko) 1/2 fi(z, k Q )A(k ) 1/2 - O^fa,) 1 / 2 , 

j = 1,2 (4.93) 

then yield 

||G,(z,fc ,±C)|| Cm x m = 0(|z|), Cer(z), j = 1,2, (4.94) 

2: — >oc 

alo ng th e ray arg(z) = it — e, u niformly with respect to Q G T(z). As a consequence 
of ( 4.94 ) and hypothesis ( 4.74 ) we have 

||Gi(z,fcb,±0-G 2 («,*o,±Ollc»x» , = 0(1*1"*). Cer(z) (4.95) 

Z — >OC 



along the ray arg(z) = 7r — e, uniformly in £ G r(z). Using ( 4.88 ), ( 4.94 ), and ( 4.95 ) 
one concludes that 



\\M lt+ (z,k )-M 2r+ (z,ko)\\cmxm = 0{\z\-^ N+t >) 



(4.96) 
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and hence by (1.76) and (4.74) that 

||Mi,_(z,to)-M 2 ,-(2,*o)||c" 



\z\ — >00 



0{\z\ 



-N\ 



(4.97) 



along the ray arg(z) = 7r — s. This proves ( 4.7q) by T heore m 4.4 
Case (iii). One observes that Mj +(z, fco) satisfy (c.f. (4.67)) 



Mj,+(z, ho) + Mj,+(z, ko) gj (z, ko)Mj,+(z, k ) = A(ko)gj(z, k a + l)A(k ), j = 1, 2. 

(4.98) 



Taking i nto a ccoun t the asymptotic e xpan sions ( 4.7£ ) a nd (4.80 ), hypothesis ( 4.75 ) 



implies ( 4.96 ) and ( 4.97 ) by Corollary 3J3. This proves ( 4.76 ) applying Theorem 4.4 
again. □ 

Without going into further details, we note that the methods employed in this 
section can easily be adapted to the case of uniformly bounded operator coefficients 
for Jacobi operators (i.e., for some C > and all k G Z, ||y4.(fc)|| B / H ) + ||-B(fc)||e(w) < 
C, A^ky 1 £ B(H), < A{k), B{k) = B(k)* , for some complex separable Hilbert 
space Ti). 

Acknowledgements. We are indebted to Steve Clark, Don Hinton, Boris Levitan, 
Mark Malamud, Fedor Rofe-Beketov, Alexander Sakhnovich, Lev Sakhnovich, and 
Barry Simon for discussions and valuable hints regarding the literature. F.G. and 
A.K. thank C. Peck and T. Tombrcllo for the hospitality of Caltech during the final 
stages of this work. 



References 

[1] Z. S. Agranovich and V.A. Marchenko, The Inverse Problem of Scattering Theory, Gordon 

and Breach, New York, 1963. 
[2] B. P. Allakhverdiev, Extensions of symmetric Schrddinger operators with matrix coefficients, 

Izvest. Math. 59, 45-62 (1995). 
[3] D. Alpay and I. Gohberg, Inverse spectral problems for difference operators with rational 

scattering matrix function, Integr. Equat. Oper. Th. 20, 125—170 (1994). 
[4] D. Alpay and I. Gohberg, Inverse spectral problem for differential operators with rational 

scattering matrix functions, J. DifT. Eq. 118, 1-19 (1995). 
[5] D. Alpay and I. Gohberg, Inverse problem for Sturm- Liouville operators with rational reflec- 
tion coefficient, Integr. Equat. Oper. Th. 30, 317-325 (1998). 
[6] A. I. Aptekarev and E. M. Nikishin, The scattering problem for a discrete Sturm- Liouville 

problem, Math. USSR Sborn. 49, 325-355 (1984). 
[7] D. Z. Arov and H. Dym, J-inner matrix functions, interpolation and inverse problems for 

canonical systems, I: foundations, Integr. equ. oper. theory 29, 373-454 (1997). 
[8] D. Z. Arov and H. Dym, J-inner matrix functions, interpolation and inverse problems for 

canonical systems, II: the inverse monodromy problem, Integr. equ. oper. theory 36, 11-70 

(2000). 

[9] D. Z. Arov and H. Dym, J-inner matrix functions, interpolation and inverse problems for 
canonical systems, III: more on the inverse monodromy problem, Integr. equ. oper. theory 
36, 127-181 (2000). 

[10] N. Asano and Y. Kato, The Marchenko equation of the Hilbert- Schmidt type for the N X N 
Zakharov-Shabat operator, Progress Theoret. Phys. 80, 774-792 (1988). 

[11] N. Asano and Y. Kato, Algebraic and Spectral Methods for Nonlinear Wave Equations, Long- 
man, New York, 1990. 

[12] F. V. Atkinson, Discrete and Continuous Boundary Problems, Academic Press, New York, 
1964. 

[13] F. V. Atkinson, On the order of magnitude of Titchmarsh- Weyl functions, Difl". Integral Eqs. 
1, 79-96 (1988). 



34 GESZTESY, KISELEV, AND MAKAROV 

[14] F. V. Atkinson, Asymptotics of the Titchmarsh- Weyl function in the matrix case, unpub- 
lished manuscript. 

[15] R. Benguria and M. Loss, A simple proof of a theorem of Laptev and Weidl, Math. Res. Lett. 
7, 195-203 (2000). 

[16] C. Bcnncwitz, A proof of the local Borg-Marchenko theorem, Comm. Math. Phys., to appear. 
[17] Ju. M. Berezanskii, On the uniqueness of the determination of Schrddinger's equation from 

its spectral function, Sov. Math. Dokl. 93, 591-594 (1953). (Russian). 
[18] Ju. M. Berezanskii, The uniqueness theorem in the inverse problem of spectral analysis for 

the Schrddinger operator, Amer. Math. Soc. Transl., Scr. 2, 35, 167-235 (1964). 
[19] Ju. Berezanskii, Expansions in Eigenf unctions of Selfadjoint Operators, Transl. Math. Mon- 

graphs, Vol. 17, Amer. Math. Soc., Providence, R.I., 1968. 
[20] Yu. M. Berezanskii, M. I. Gekhtman, and M. E. Shmoish, Integration of some chains of 

nonlinear difference equations by the method of the inverse spectral problem, Ukrain. Math. 

J. 38, 74-78 (1986). 

[21] Yu. M. Berezanskii and M. I. Gekhtman, Inverse problem of the spectral analysis and non- 

abelian chains of nonlinear equations, Ukrain. Math. J. 42, 645—658 (1990). 
[22] R. Bhatia, Matrix Analysis, Springer, New York, 1997. 

[23] R. Bhatia, C. Davies, and A. Mcintosh, Perturbation of spectral subspaces and solution of 
linear operator equations, Linear Algebra Appls. 52/53, 45-67 (1983). 

[24] G. Borg, Uniqueness theorems in the spectral theory of y" + (X—q(x))y = 0, Proc. 11th Scandi- 
navian Congress of Mathematicians, Johan Grundt Tanums Forlag, Oslo, 1952, pp. 276-287. 

[25] R. Carlson, Large eigenvalues and trace formulas for matrix Sturm- Liouville problems, SIAM 
J. Math. Anal. 30, 949-962 (1999). 

[26] R. Carlson, Compactness of Floquet isospectral sets for the matrix Hill's equation, Proc. 
Amer. Math. Soc. 128, 2933-2941 (2000). 

[27] R. Carlson, Eigenvalue estimates and trace formulas for the matrix Hill's equation, J. Diff. 
Eq. 167, 211-244 (2000). 

[28] R. Carlson, An inverse problem for the matrix Schrddinger equation, preprint, 2000. 

[29] R. Carmona and J. Lacroix, Spectral Theory of Random Schrddinger operators, Birkhauscr, 
Boston, MA, 1990. 

[30] I. Chcrcdnik, Basic Methods of Soliton Theory, World Scientific, Singapore, 1996. 
[31] H.-H. Chern, On the construction of isospectral vectorial Sturm- Liouville differential equa- 
tions, preprint, 1998. 

[32] H.-H. Chern and C-L. Shen, On the n-dimensional Ambarzumyan's theorem, Inverse Problems 
13, 15-18 (1997). 

[33] S. L. Clark, A criterion for absolute continuity of the continuous spectrum of a Hamiltonian 

system, J. Math. Anal. Appl. 151, 108-128 (1990). 
[34] S. L. Clark, Asymptotic behavior of the Titchmarsh- Weyl coefficient for a coupled second 

order system, in ^Ordinary and Delay Differential Equations" , J. Wiener and J. K. Hale 

(cds.), Longman, New York, 1992, pp. 24-28. 
[35] S. L. Clark, On the absolutely continuous spectrum of a vector-matrix Dirac system, Proc. 

Roy. Soc. Edinburgh 124A, 253-262 (1994). 
[36] S. Clark and F. Gesztesy, Weyl- Titchmarsh M -function asymptotics for matrix-valued 

Schrddinger operators, Proc. London Math. Soc, to appear. 
[37] S. Clark and F. Gesztesy, Weyl- Titchmarsh M -function asymptotics, local uniqueness results, 

trace formulas, and Borg-type theorems for Dirac operators, preprint, 2001. 
[38] S. Clark, F. Gesztesy, H. Holden, and B. M. Levitan, Borg-type theorems for matrix-valued 

Schrddinger and Dirac operators, J. Diff. Eqs. 167, 181-210 (2000). 
[39] S. Clark and D. Hinton, A Liapunov inequality for linear Hamiltonian systems, Math. Inequ. 

Appl. 1, 201-209 (1998). 
[40] P. Delsarte, Y. V. Genin, and Y. G. Kamp, Orthogonal polynomial matrices on the unit 

circle, IEEE Trans. Circ. Syst. 25, 149-160 (1978). 
[41] V. I. Dcrguzov, The spectrum of Hamilton's operator with periodic coefficients, Vcstnik 

Leningrad Univ. Math. 12, 280-285 (1980). 
[42] B. Despres, The Borg theorem for the vectorial Hill's equation, Inverse Probl. 11, 97-121 

(1995). 

[43] L. A. Dickey, Soliton Equations and Hamiltonian Systems, World Scientific, Singapore, 1991. 



UNIQUENESS FOR SCHRODINGER, JACOBI, AND DIRAC-TYPE OPERATORS 



35 



[44] B. A. Dubrovin, Completely integrable Hamiltonian Systems associated with Matrix operators 

and Abelian varieties, Funct. Anal. Appl. 11, 265-277 (1977). 
[45] B. A. Dubrovin, Matrix finite- zone operators, Revs. Sci. Tech. 23, 20-50 (1983). 
[46] A. J. Duran and W. Van Assche, Orthogonal matrix polynomials and higher-order recurrence 

relations, Lin. Algebra Appl. 219, 261-280 (1995). 
[47] A. J. Duran and P. Lopez-Rodriguez, Orthogonal matrix polynomials: zeroe and Blumenthal's 

theorem, J. Approx. Th. 84, 96-118 (1996). 
[48] M. Fukushima, A spectral representation on ordinary linear difference equation with operator- 
valued coefficients of the second order, J. Math. Phys. 17, 1084-1072 (1976). 
[49] M. G. Gasymov, The inverse scattering problem for a system of Dirac equations of order 2n, 

Trans. Moscow Math. Soc. 19, 41-119 (1968). 
[50] M. G. Gasymov and B. M. Levitan, The inverse problem for a Dirac system, Sov. Math. 

Dokl. 7, 495-499 (1966). 
[51] I. M. Gcl'fand and L. A. Dikii, The resolvent and Hamiltonian systems, Funct. Anal. Appl. 

11, 93-105 (1977). 

[52] I. M. GeFfand and V. B. Lidskij, On the structure of the regions of stability of linear canonical 
systems of differential equations with periodic coefficients, in "Izrail M. Gelfand, Collected 
Papers", Vol. I, S. G. Gindikin, V. W. Guillemin, A. A. Kirillov, B. Kostant, S. Sternberg 
(cds.), Springer, Berlin, 1987, pp. 466-504. (Amcr. Math. Soc. Transl. (2) 8, 143-181 (1958).) 

[53] J. S. Gcronimo, Matrix orthogonal polynomials on the unit circle, J. Math. Phys. 22, 1359- 
1365 (1981). 

[54] J. S. Gcronimo, Scattering theory and matrix orthogonal polynomials on the real line, Circuits 
Syst. Signal Process. 1, 471-495 (1982). 

[55] F. Gesztesy and H. Holdcn, On trace formulas for Schrddinger-type operators, in Multi- 
particle Quantum Scattering with Applications to Nuclear, Atomic and Molecular Physics, 
D. G. Truhlar and B. Simon (eds.), Springer, New York, 1997, pp. 121-145. 

[56] F. Gesztesy, K. A. Makarov, and S. N. Naboko, The spectral shift operator, in Mathematical 
Results in Quantum Mechanics, J. Dittrich, P. Exner, and M. Tater (eds.), Operator Theory: 
Advances and Applications, Vol. 108, Birkhauser, Basel, 1999, pp. 59-90. 

[57] F. Gesztesy and B. Simon, Uniqueness theorems in inverse spectral theory for one-dimensinal 
Schrddinger operators, Trans. Amer. Math. Soc. 348, 349-373 (1996). 

[58] F. Gesztesy and B. Simon, The § function, Acta Math. 176, 49-71 (1996). 

[59] F. Gesztesy and B. Simon, A new approach to inverse spectral theory, II. General real po- 
tentials and the connection to the spectral measure, Ann. of Math. 152, 593—643 (2000). 

[60] F. Gesztesy and B. Simon, On local Borg-Marchenko uniqueness results, Commun. Math. 
Phys. 211, 273-287 (2000). 

[61] F. Gesztesy and E. Tsckanovskii, On matrix-valued Herglotz functions, Math. Nachr. 218, 
61-138 (2000). 

[62] I. Gohbcrg, M. A. Kaashoek, and A. L. Sakhnovich, Sturm- Liouville systems with rational 
Weyl functions: explicit formulas and applications, Integral Eq. Operator Th. 30, 338-377, 
(1998). 

[63] K. E. Gustafson and D. K. M. Rao, Numerical Range: The Field of Values of Linear Oper- 
ators and Matrices, Springer, New York, 1997. 

[64] E. Gutkin and R. Johnson, Intersection theory for linear eigenvalue problems, J. reinc angew. 
Math. 401, 1-24 (1989). 

[65] E. Hille, Lectures on Ordinary Differential Equations, Addison-Wesley, Reading, 1969. 

[66] D. B. Hinton and A. Schneider, On the Titchmarsh-Weyl coefficients for singular S- 
Hermitian Systems I, Math. Nachr. 163, 323-342 (1993). 

[67] D. B. Hinton and A. Schneider, On the Titchmarsh-Weyl coefficients for singular S- 
Hermitian Systems II, Math. Nachr. 185, 67-84 (1997). 

[68] D. B. Hinton and J. K. Shaw, On Titchmarsh-Weyl M(X)- functions for linear Hamiltonian 
systems, J. Diff. Eqs. 40, 316-342 (1981). 

[69] D. B. Hinton and J. K. Shaw, On the spectrum of a singular Hamiltonian system, Quacst. 
Math. 5, 29-81 (1982). 

[70] D. B. Hinton and J. K. Shaw, Hamiltonian systems of limit point or limit circle type with 

both endpoints singular, J. Diff. Eqs. 50, 444-464 (1983). 
[71] D. B. Hinton and J. K. Shaw, On boundary value problems for Hamiltonian systems with 

two singular points, SIAM J. Math. Anal. 15, 272-286 (1984). 



36 GESZTESY, KISELEV, AND MAKAROV 

[72] D. B. Hinton and J. K. Shaw, On the spectrum of a singular Hamiltonian system, II, Quacst. 
Math. 10, 1-48 (1986). 

[73] M. Jodeit and B. M. Levitan, Isospectral vector-valued Sturm- Liouville problems, Lett. Math. 

Phys. 43, 117-122 (1998). 
[74] M. Jodeit and B. M. Levitan, The isospectrality problem for some vector-valued Sturm- 

Liouville boundary problems, preprint, 1999. 
[75] R. A. Johnson, m-Functions and Floquet exponents for linear differential systems, Ann. Mat. 

Pura Appl., Ser. 4, 147, 211-248 (1987). 
[76] R. Johnson and J. Moser, The rotation number for almost periodic potentials, Commun. 

Math. Phys. 84, 403-438 (1982). 
[77] R. Johnson, S. Novo, and R. Obaya, Ergodic properties and Weyl M -functions for random 

linear Hamiltonian systems, Proc. Roy. Soc. Edinburgh 130A, 1045-1079 (2000). 
[78] V. I. Kogan and F. S. Rofc-Bckctov, On square-integrable solutions of symmetric systems of 

differential equations of arbitrary order, Proc. Roy. Soc. Edinburgh 74A, 1—40 (1974). 
[79] S. Kotani, Ljapunov indices determine absolutely continuous spectra of stationary random 

one- dimensional Schrddinger operators, in Stochastic Analysis, K. Ito (ed.), North-Holland, 

Amsterdam, 1984, pp. 225-247. 
[80] S. Kotani, On an inverse problem for random Schrddinger operators, Contemporary Math. 

41, 267-281 (1985). 

[81] S. Kotani, One- dimensional random Schrddinger operators and Herglotz functions, in Prob- 
abilistic Methods in Mathematical Physics, K. Ito and N. Ikeda (eds.), Academic Press, New 
York, 1987, pp. 219-250. 

[82] S. Kotani and B. Simon, Stochastic Schrddinger operators and Jacobi matrices on the strip, 
Commun. Math. Phys. 119, 403-429 (1988). 

[83] A. M. Krall, M(A) theory for singular Hamiltonian systems with one singular point, SIAM 
J. Math. Anal. 20, 664-700 (1989). 

[84] A. M. Krall, M{\) theory for singular Hamiltonian systems with two singular points, SIAM 
J. Math. Anal. 20, 701-715 (1989). 

[85] M. G. Krein, Foundations of the theory of \-zones of stability of a canonical system of 
linear differential equations with periodic coefficients, Amer. Math. Soc. Transl. (2) 120, 
1-70 (1983). 

[86] M. G. Krein, On tests for stable boundedness of solutions of periodic canonical systems, 

Amer. Math. Soc. Transl. (2) 120, 71-110 (1983). 
[87] J. Lacroix, The random Schrddinger operator in a strip, in Probability Measures on Groups 

VII, H. Heyer (ed.), Lecture Notes in Math. 1064, Springer, Berlin, 1984, pp. 280-297. 
[88] J. Lacroix, Computations of the sum of positive Lyapunov exponents for the Llyod model in 

a strip, in Lyapunov Exponents, L. Arnold and V. Wihstutz (eds.), Lecture Notes in Math. 

1186, Springer, Berlin, 1986, pp. 258-264. 
[89] A. Laptev and T. Weidl, Sharp Lieb-Thirring inequalities in high dimensions, Acta Math. 

184, 87-111 (2000). 

[90] Z. L. Leibenzon, The uniqueness of the solution of the inverse problem for ordinary differ- 
ential operators of order n > 2 and the transformation of such operators, Sov. Math. Dokl. 
3, 100-104 (1962). 

[91] Z. L. Leibenzon, A connection between the inverse problem and the completeness of eigen- 

functions, Sov. Math. Dokl. 3, 1045-1048 (1962). 
[92] M. Lesch and M. Malamud, The inverse spectral problem for first order systems on the 

half line, Operator Theory: Advances and Applications, Vol. 117, Birkhauser, Basel, 2000, 

pp. 199-238. 

[93] M. Lesch and M. Malamud, On the number of square integrable solutions and self-adjointness 
of symmetric first order systems of differential equations, preprint, 2000. 

[94] B. M. Levitan and I. S. Sargsjan, Introduction to Spectral Theory, Amer. Math. Soc, Provi- 
dence, R. I., 1975. 

[95] B. M. Levitan and M. Otelbaev, Conditions for self-adjointness of the Schrddinger and Dirac 

operators, Trans. Moscow. Math. Soc. 42, 139-156 (1982). 
[96] B. M. Levitan and I. S. Sargsjan, Sturm- Liouville and Dirac Operators, Kluwer, Dordrecht, 

1991. 

[97] P. Lopez-Rodriguez, Riesz's theorem for orthogonal matrix polynomials, Constr. Approx. 15, 
135-151 (1999). 



UNIQUENESS FOR SCHRODINGER, JACOBI, AND DIRAC-TYPE OPERATORS 



37 



[98] F. G. Maksudov, E. M. Bairamov, and R. U. Orudzheva, The inverse scattering problem for 
an infinite Jacobi matrix with operator elements, Russ. Acad. Sci. Dokl. Math. 45, 366-370 
(1992). 

[99] M. M. Malamud, Similarity of Volterra operators and related questions of the theory of 
differential equations of fractional order, Trans. Moscow Math. Soc. 55, 57—122 (1994). 

[100] M. M. Malamud, A connection between the potential matrix of the Dirac system and its 
Wronskian, Dokl. Math. 52, 296-299 (1995). 

[101] M. M. Malamud, Uniqueness questions in inverse problems for systems of ordinary differ- 
ential equations on a finite interval, Trans. Moscow Math. Soc. 60, 173-224 (1999). 

[102] M. M. Malamud, Borg type theorems for first-order systems on a finite interval, Funct. 
Anal. Appl. 33, 64-68 (1999). 

[103] Yu. I. Manin, Matrix solitons and bundles over curves with singularities, Funct. Anal. Appl. 
12, 286-295 (1978). 

[104] V. A. Marchcnko, Certain problems in the theory of second-order differential operators, 

Doklady Akad. Nauk SSSR 72, 457-460 (1950) (Russian). 
[105] V. A. Marcenko, Some questions in the theory of one- dimensional linear differential opera- 
tors of the second order. I, Trudy Moskov. Mat. Obsc. 1, 327-420 (1952) (Russian); English 

transl. in Amcr. Math. Soc. Transl. (2) 101, 1-104 (1973). 
[106] V. A. Marchcnko, Sturm- Liouville Operators and Applications, Birkhauser, Basel, 1986. 
[107] V. A. Marchcnko, Nonlinear Equations and Operator Algebras, Reidel, Dordrecht, 1988. 
[108] L. Martinez Alonso and E. Olmedilla, Trace identities in the inverse scattering transform 

method associated with matrix Schrddinger operators, J. Math. Phys. 23, 2116-2121 (1982). 
[109] V. V. Martynov, Condition for discreteness and continuity of the spectrum of a selfadjoint 

system of first- order differential equations, Sov. Math. Dokl. 6, 1546-1550 (1965). 
[110] R. Mcnnickcn, A. L. Sakhnovich, and C. Tretter, Direct and inverse spectral problem for a 

system of differential equations depending rationally on the spectral parameter, Duke Math. 

J., to appear. 

[Ill] R. G. Newton and R. Jost, The construction of potentials from the S -matrix for systems of 

differential equations, Nuovo Cim. 1, 590-622 (1955). 
[112] E. Olmedilla, L. Martinez Alonso, and F. Guil, Infinite-dimensional Hamiltonian systems 

associated with matrix Schrddinger operators, Nuovo Cim. 61 B, 49-61 (1981). 
[113] S. A. Orlov, Nested matrix disks analytically depending on a parameter, and theorems on 

the invariance of ranks of radii of limiting disks, Math. USSR Izv. 10, 565-613 (1976). 
[114] V. G. Papanicolaou, Trace formulas and the behavior of large eigenvalues, SIAM J. Math. 

Anal. 26, 218-237 (1995). 
[115] V. Q. Phong, The operator equation AX — XB = C with unbounded operators A and B 

and related abstract Cauchy problems, Math. Z. 208, 567-588 (1991). 
[116] F. S. Rofe-Beketov, Expansions in eigenf unctions of infinite systems of differential equations 

in the non-self-adjoint and self-adjoint cases, Mat. Sb. 51, 293-342 (1960). (Russian). 
[117] F. S. Rofe-Bckctov, The spectrum of non-selfadjoint differential operators with periodic 

coefficients, Sov. Math. Dokl. 4, 1563-1566 (1963). 
[118] F. S. Rofe-Beketov, Selfadjoint extensions of differential operators in a space of vector 

functions, Sov. Math. Dokl. 10, 188-192 (1969). 
[119] F. S. Rofe-Bckctov and A. M. HoFkin, On the connection between spectral and oscillation 

properties of the Sturm- Liouville matrix problem, Math. USSR Sbornik 31, 365-378 (1977). 
[120] A. L. Sakhnovich, Asymptotic behavior of spectral functions of an S-node, Sov. Math. (Iz. 

VUZ) 32:9, 92-105 (1988). 
[121] A. L. Sakhnovich, Nonlinear Schrddinger equation on a semi-axis and an inverse problem 

associated with it, Ukrain. Math. J. 42, 316-323 (1990). 
[122] A. L. Sakhnovich, Spectral functions of a canonical system of order 2n, Math. USSR Sbornik 

71, 355-369 (1992). 

[123] A. Sakhnovich, Iterated Backlund-Darboux transform for canonical systems, J. Funct. Anal. 
144, 359-370 (1997). 

[124] A. Sakhnovich, Canonical systems and transfer matrix-functions, Proc. Amer. Math. Soc. 
125, 1451-1455 (1997). 

[125] A. Sakhnovich, Dirac type and canonical systems: spectral functions, direct and inverse 
problems, preprint, 2001. 



:!8 



GESZTESY, KISELEV, AND MAKAROV 



[126; 
[127; 

[128; 
[129 

[130; 

[131 

[132; 
[133; 
[134; 

[135 
[136 
[137; 

[138; 
[139; 
[140; 

[141 
[142; 
[143; 
[144 
[145; 



L. A. Sakhnovich, Evolution of spectral data and nonlinear equations, Ukrain. Math. J. 40, 
459-461 (1988). 

L. A. Sakhnovich, Inverse problems for equations systems, in Matrix and Operator 
Valued Functions: The Vladimir Petrovich Potapov Memorial Volume, I. Gohberg and 
L. A. Sakhnovich (eds.), Operator Theory: Advances and Applications, Vol. 72, Birkhauser, 
Basel, 1994, pp. 202-211. 

L. A. Sakhnovich, Method of operator identities and problems of analysis, St. Petersburg 
Math. J. 5, 1-69 (1994). 

L. A. Sakhnovich, Spectral problem on half-axis, Methods Funct. Anal. Topology 2, 128-140 
(1996). 

L. A. Sakhnovich, Interpolation Theory and its Applications, Kluwer, Dordrecht, 1997. 

L. A. Sakhnovich, Spectral analysis of a class of canonical differential systems, St. Peters- 
burg Math. J. 10, 147-158 (1999). 

L. A. Sakhnovich, Spectral Theory of Canonical Differential Systems. Method of Operator 
Identities, Operator Theory: Advances and Applications, Vol. 107, Birkhauser, Basel, 1999. 

C.-L. Shen and C.-T. Shieh, Two inverse eigenvalue problems for vectorial Sturm- Liouville 
equations, Inverse Probl. 14, 1331-1343 (1998). 

B. Simon, A new aproach to inverse spectral theory, I. Fundamental formalism, Ann. of 
Math. 150, 1029-1057 (1999). 

G. Teschl, Trace formulas and inverse spectral theory for Jacobi operators, Commun. Math. 
Phys. 196, 175-202 (1998). 

G. Teschl, Jacobi Operators and Completely Integrable Nonlinear Lattices, Mathematical 
Surveys and Monographs, Vol. 72, Amer. Math. Soc, Providence, RI, 2000. 

I. Trooshin, Asymptotics for the spectral and Weyl functions of the operator-valued Sturm- 
Liouville problem, in "Inverse Problems and Related Topics", G. Nakamura, S. Saitoh, 
J. K. Seo, and M. Yamamoto (eds.), Chapman & Hall/CRC, Res. Notes Math. 419, Boca 
Raton, FL, 2000, pp. 189-208. 

M. Wadati and T. Kamijo, On the extension of inverse scattering method, Progr. Theoret. 
Phys. 52, 397-414 (1974). 

J. Weidmann, Oszillationsmethoden fur Systeme gewdhnlicher Differentialgleichungen, 
Math. Z. 119, 349-373 (1971). 

J. Weidmann, Spectral Theory of Ordinary Differential Operators, Lecture Notes in Math. 
1258, Springer, Berlin, 1987. 

W. Weyl, Uber gewdhnliche Differentialgleichungen mit Singularitaten und die zugehdren 
Entwicklungen willkurlicher Funktionen, Math. Ann. 68, 220—269 (1910). 

V. A. Yakubovich, Nonoscillation of linear periodic Hamiltonian equations, and related 
topics, St. Petersburg Math. J. 3, 1165-1188 (1992). 

V. A. Yakubovich and V. M. Starzhinskii, Linear Differential Equations with Periodic Co- 
efficients 1, Wiley, New York, 1975. 

V. A. Yakubovich and V. M. Starzhinskii, Linear Differential Equations with Periodic Co- 
efficients 2, Wiley, New York, 1975. 

X. Zhou, Inverse scattering transform for systems with rational spectral dependence, J. Diff. 
Eqs. 115, 277-303 (1995). 



Department of Mathematics, University of Missouri, Columbia, MO 65211, USA 
E-mai address: fritzOmath.missouri.edu 



URL: http: / /www. math. mis sour i . edu/people/f gesztesy .html 



Department of Mathematics, University of Missouri, Columbia, MO 65211, USA 
E-mai address: makarov@math.missouri.edu 



URL: http: / /www. math. mis sour i . edu/people/kmakarov. html 



Department of Mathematics, University of Chicago, 5734 South University Avenue, 
Chicago, IL 60637-1546 USA 

E-mail address: kiselev9math.uchicago.edu 



